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High resolution optical spectra of HD 37903 and HD 147009, which illuminate 

the reflection nebulae, NGC 2023 and vdB 102, were obtained for comparison 
with our results for HD 200775 and NGC 7023. Ground-based measurements of 
the molecules, CH, C2, and CN, and the atoms, Na I and K I, were analyzed 
to extract physical conditions in the foreground cloud. Estimates of the gas 
density, gas temperature and flux of ultraviolet radiation were derived and were 
compared with the results from infrared and radio studies of the main molecular 
cloud. The conditions are similar to those found in studies of diffuse clouds. 
The foreground material is less dense than the gas in the molecular cloud behind 
the star(s). The gas temperature was set at 40 K, the temperature determined 
for the foreground gas in NGC 7023. The flux of ultraviolet radiation was 
found to be less intense than in the molecular material behind the star(s). 
The column densities of Na I and K I were reproduced reasonably well when 
the extinction curve for the specific line of sight was adopted. We obtained 
NEWSIPS data from the lUE archive for HD 37903 and HD 200775. The 
ultraviolet data on C I and CO allow extraction of the physical conditions by 
alternate methods. General agreement among the various diagnostics was found, 
leading to self-consistent pictures of the foreground photodissociation regions. 

An Appendix describes checks on the usefulness of lUE NEWSIPS data for 
interstellar studies. Equivalent widths are compared with a previous analysis of 
lUE observations for interstellar gas toward 20 Aquilae. Excellent agreement is 
found with NEWSIPS results having smaller errors. A comparison of NEWSIPS 
data for C I toward X Per with data acquired with the Hubble Space Telescope 
shows similar levels of agreement in equivalent widths and that the derived 
column densities from lU E results are accurate to better than a factor of 2 for 
absorption lines of moderate strength. 

Subject headings: ISM: atoms — molecules — refiection nebulae — ISM: 
individual (NGC 2023, vdB 102, and NGC 7023) — Stars: individual 
(HD 37903, HD 147009, 20 Aquilae, and HD 200775) 

1. Introduction 

Stars form under the gravitational collapse of gas and dust in interstellar molecular 
clouds. When a newly-formed star interacts with the surrounding material, the cloud's 
physical conditions are altered. For instance, early type stars create photodissociation 



-3- 



regions (PDRs), which are boundary layers between a source of far ultraviolet (FUV) 
radiation and the predominantly neutral cloud (Jaffe et al. 1990). These boundaries also 
exist between H II regions produced by O stars and the ambient material surrounding 
them. In photodissociation regions, the FUV radiation dominates the heating and chemical 
processes (Tielens & Hollenbach 1985a, 1985b; Hollenbach & Tielens 1999, HT99) and 
controls which atoms and molecules are observed. The observed relative abundances of 
atomic and molecular species and the use of chemical models yields the physical conditions 
in the PDR. These conditions provide clues to processes that occurred immediately following 
star formation. In this paper, we describe results for foreground PDR's in three reflection 
nebulae based on measurements of absorption at visible and ultraviolet wavelengths 
occurring at the velocity of the PDR behind the illuminating source seen via emission 
lines. While absorption is sometimes seen at other velocities, the main components all have 
velocities consistent with the molecular cloud. From this fact, we infer an association with 
the material behind the star. 

PDRs in molecular clouds and reflection nebulae have been studied extensively 
at all wavelengths. Analyses have been based on measurements from x-rays to radio 
wavelengths. Speciflc observations on the nebulae studied here are given below. In an effort 
to obtain a comprehensive understanding of the conditions that prevail in these regions, 
one dimensional, time-independent models have been developed (e.g., Tielens & Hollenbach 
1985a, 1985b; van Dishoeck & Black 1986; Hollenbach, Takahashi & Tielens 1991; Draine & 
Bertoldi 1996, 1999). More reahstic PDR models with clumpy structure and time dependent 
effects have been attempted as well (Meixner & Tielens 1993; Hollenbach & Natta 1995; 
Bertoldi & Draine 1996). Hollenbach & Tielens (1999) provide a comprehensive review of 
observations and modeling efforts for photodissociation regions. 

Little attention has been given to studying PDRs through atomic and molecular 
absorption in the visible and ultraviolet portions of the spectrum observed along the line 
of sight to the central or illuminating star(s) . Federman et al. (1997 - hereafter Paper I) 
analyzed visible spectra for the line of sight toward HD 200775, which illuminates the 
reflection nebula NGC 7023. In an effort to expand upon this work, we obtained high 
resolution visible spectra toward NGC 2023 and vdB 102 as well as ultraviolet data on atoms 
and CO molecules in the reflection nebulae, NGC 2023 and NGC 7023. The ultraviolet 
data were obtained from the International Ultraviolet Explorer {lUE) archive, which was 
recalibrated with the New Spectral Image Processing System (NEWSIPS). The visible and 
UV absorption spectra yield information about the chemical composition and the physical 
conditions present in the foreground PDR. A more complete picture of PDRs emerges when 
the physical conditions in front of the star are compared to those inferred from analyses of 
data obtained at other wavelengths of the main molecular cloud behind the star. 
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1.1. NGC 2023 

NGC 2023 is a reflection nebula embedded in the L1630 molecular cloud (Jaffe et al. 
1990) in the constellation of Orion. Infrared images of the nebula (Gatley et al. 1987; 
Field et al. 1994, 1998) show that it has a 6' extent with visible filamentary structures 
throughout and is approximately centered on the star HD 37903. In a far- infrared study, 
Harvey, Thronson, & Gatley (1980) showed that the molecular cloud lies behind HD 37903, 
which is located at a distance of 470 pc (Ferryman et al. 1997). Several Herbig-Haro 
objects, which arc associated with pre-main sequence stars, are observed in the nebula. Two 
Herbig-Haro objects (HH4 and HH5) in the southeastern part of the nebula are illuminated 
by star C (Malin, Ogura & Walsh 1987). Star C is thought to be a T Tauri star based 
on spectroscopy and its near infrared excess. Many embedded infrared sources have also 
been detected (Strom et al. 1975). Wyrowski et al. (2000) detected a ultra-high density 
clump {n ~ 10^ cm~^) which they propose could be the early stages of a protostar. Thus, 
NGC 2023 is an active region of star formation (Strom et al. 1975; Mahn et al. 1987; Lada 
et al. 1991; Preyberg & Schmitt 1995). 

The primary illuminating source for this nebula is HD 37903, a B1.5 V star with an 
effective temperature of approximately 22,000 K (Field et al. 1998). The radiation from 
this star has enough energy at A < 912 A to form an H II region, 0.015 pc in size, detected 
by radio continuum emission at the star's position within the nebula (Pankonin & Walmsley 
1976, 1978; Wyrowski et al. 2000). The photons with A > 912 A create a photodissociation 
region between the H II region and the rest of the neutral cloud. The proximity of the 
nebula and its bubble shaped geometry (Gatley et al. 1987) make it an ideal target for 
studying the interaction of the UV radiation from HD 37903 with the molecular material 
in L1630. tHE refiection nebula and its PDR have been extensively studied via x-ray 
(Preyberg & Schmitt 1995), ultraviolet (Witt, Bohhn, & Stecher 1984 - hereafter WBS84; 
CardeUi, Clayton, & Mathis 1989 - hereafter CCM89; Buss et al. 1994), infrared (e.g.. 
Martini, Sellgren, & DePoy 1999 and references therein), millimeter (Jaffe et al. 1990), 
and radio observations (Puente, Martm-Pintado & Gaume 1995; Lebron & Rodriguez 1997; 
Wyrowski et al. 1997, 2000). The studies relevant to our work will be discussed in detail 
later. 



1.2. vdB 102 

The central star of vdB 102 is HD 147009, an AO V star, with an effective temperature 
of about 9,500 K (de Geus, de Zeeuw & Lub 1989 - hereafter GZL89). This nebula was 
first catalogued and photometrically studied (van den Bergh 1966; Racine 1968) in the late 
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1960's. The nebula is 160 pc away (Perryman et al. 1997) and is one of 10 reflection nebulae 
in the Sco 0B2 Association (van den Bergh 1966; Vallee 1987; Eggen 1998). The stars in 
Sco 0B2 illuminate different portions of the same molecular cloud. Radio observations 
reveal the presence of the parent molecular cloud (Kutner et al. 1980; Vallee 1987). Few 
observational efforts have examined the reflection nebula and its illuminating star (Chaffee 
& White 1982; Cappa de Nicolau & P5ppcl 1986; Vallee 1987; Federman et al. 1994). 
For this reason, the present study greatly improves upon our knowledge of PDRs around 
relatively cool stars. 

The remainder of the paper is organized in the following manner. In §2 we discuss 
the optical and ultraviolet observations obtained for this study, as well as the consistency 
of our measurements with previous efforts. We describe the chemical and atomic analyses 
performed here and the results of these analyses as well as a comparison with our earlier 
results on HD 200775 (Paper I) in §3. We discuss in §4 the derived physical conditions 
and compare them with the physical conditions obtained through diagnostics at other 
wavelengths. Our conclusions are found in §5. An Appendix describes checks on the 
usefulness of NEWSIPS data for interstellar studies such as ours. 



2. Observations and Data Reduction 
2.1. Ground-based Observations 

2.1.1. McDonald 2dcoude Data 

High-resolution echeUe spectra of HD 37903, HD 147009 and 20 Aql were taken 
with the 2.7 m telescope of the University of Texas McDonald Observatory during June 
1996, July 1997, and January 1998. Observations of 20 Aql were made since it allows 

for the comparison of our optical data with others and it is located behind a relatively 
isolated cloud for an additional test of the atomic analysis discussed below. The "2dcoude" 
spectrograph (Tull et al. 1995) was used with 2 camera setups and a Tektronix charge 
coupled device (CCD) was the detector. One camera provided high spectral resolution with 
a resolving power at 4000 A of R ~ 180,000 (Av ~ 2 km s^^) and covered the wavelength 
range 3800 A to 4400 A. This high resolution mode yielded spectra containing 8 disjoint 
echelle orders approximately 15 A wide, separated by approximately 65 A. This setup 
provided spectra on the following interstellar species CH+, CH, CN, Ca II K, and Ca I. 
The second camera setting resulted in moderate spectral resolution with a resolving power 
at 7000 A of R ?5i 52,000 (Av ~ 6 km s~^) and covered a wavelength range from 3600 
A to 9700 A. The moderate resolution spectra consisted of 59 orders with an overlap of 
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about 15 A at the blue end. The inter-order separation increased with wavelength until 
there was approximately 150 A between orders at near infrared wavelengths. In addition 
to observing the same species obtained from the high resolution spectra, we also obtained 
data on Ca II H, Na I Di, K I, and C2. The resolution was determined from the full width 
at half maximum (FWHM) of Thorium- Argon (Th-Ar) lines. 

The data were reduced in a standard way utilizing the NO AO SUN/IRAF software 
(Revision 2.10.4). Dark, bias, and fiat lamp exposures were taken each night to remove any 
instrumental effects due to the CCD detector. Comparison spectra were taken periodically 
throughout the night, typically every two hours. The average bias exposure was subtracted 
from all raw stellar, comparison (Th-Ar) and fiat images. The scattered light was fitted 
by a low order polynomial, in both the dispersion direction and perpendicular to it for 
the multi-order observations, and removed. The pixel to pixel sensitivity was removed by 
dividing the normalized average flat into the stellar spectra. Next, the pixels perpendicular 
to the dispersion were summed in each order for each stellar and comparison lamp exposure. 
The extracted spectra were placed on an appropriate wavelength scale with the Th-Ar 
comparison spectra, and Doppler-corrected. The spectra were coadded and normalized to 
unity yielding a final spectrum with high signal to noise (based on the root-mean-square 
deviations in the continuum), typically 100:1. 

2.1.2. McDonald Observations with the 6-Foot Camera 

The 6-foot camera of the coude spectrograph on the 2.7 m telescope was used to 
acquire high-resolution spectra of HD 147009 and 20 Aql. Observations, at R ~ 200,000, of 
interstellar Na I D AA5896, 5890, K I A7699, CH+ A4232, and Ca I A4226 were obtained in 
June 1990, Ca II K A3933 observations were made in June 1991 and CH A4300 observations 
were obtained in June 1992. The TI2 CCD with 15 // pixels was used in all setups. Since 
no cross disperser was used, single orders of the echelle grating were detected in the vicinity 
of the line of interest. Appropriate interference filters were placed at the entrance slit to 
remove light from unwanted orders. Scattered light was found to be negligible since only 
one order was observed. Comparison spectra from a Th-Ar hollow cathode were taken 
throughout the night, while biases and fiat fields were acquired each night. Dark frames 
were also obtained as a check on the amount of thermal noise from the cooled CCD. The slit 
was imaged onto 2.5 pixels for all setups, providing a resolution of 1.6 to 1.8 km s~^. The 
data were reduced in the standard manner described above. The resulting signal-to-noise 
ratios ranged from 30 to 200. 



- 7- 



2.1.3. CTIO Observations 

The stars, HD 147009 and 20 Aql, were observed with the Bench Mounted Echelle on 
the 1.5 m telescope in July 1991. The setup involved the 590 mm long camera, the Ritchey 
Chretien grating from the 4-m telescope as a cross disperser, the Til 800 x 800 pixel CCD, 
and a CUSO4 blocking filter. The entrance sht at the end of the 200 //m fiber optic feed 
was set at 45 /xm, illuminating 2.6 pixels as determined by the full width at half maximum 
of Th-Ar lines. The resulting instrumental resolving power was about 65,000. The setup 
allowed measurement of lines between 3800 A and 4400 A. Bias frames and flat fields 
were acquired each night. The wavelength scale was determined from spectra of a Th-Ar 
hollow cathode taken throughout the night. The data were reduced in the standard manner 
described above. The extracted spectra had signal-to-noise ratios between 40 and 150. 

A sample of the rectified ground-based spectra for atoms and molecules toward 
HD 37903 are shown in Figures 1 — 3 and Figures 4 — 6 exhibit spectra of HD 147009. 
The equivalent width, Wx, for a line was determined by a Gaussian fit; single Gaussians 
were adequate for all except the atomic species Na I and Ca II toward HD 37903. The 
measured equivalent widths and vlsr for the ground-based results toward HD 37903 and 
HD 147009 are shown in Tables 1 — 3. Where applicable, measurements taken from the 
literature are also shown. All WxS for a given line agree within measured uncertainties. 
The differences in vlsr are attributed to differences in spectral resolution and to systematic 
effects involving the wavelength solution. Where no interstellar absorption was detected, a 
3-0" upper limit is given. 

Each W\ was converted into a column density utilizing curves of growth. The column 
densities for each species were consistent within the errors from individual measurements. 
The final column densities reported in Tables 4 and 5 for HD 37903 and HD 147009, 
respectively, are the weighted average of all detections reported in this work and the 
literature. For the C2 results reported in Table 3, the most restrictive upper limits for a 
rotational level were summed to get the total upper limit in column density. 

The Doppler parameter (6- value) was set at 1 km s~^ for all lines except those of Ca II 
and CH"*", where a lvalue of 2.5 km was used (e.g., Welty, Morton, & Hobbs 1996). For 
Na I and K I absorption associated with the molecular material of the FDR, 6-values of 2.0 
km s"^ toward HD 37903 and 1.5 km s"^ toward HD 147009 were required by the atomic 
analysis discussed below. These 6-values are similar to other sightlines studied in Orion and 
the Sco 0B2 association based on high resolution surveys of interstellar Na I Di and K I 
(Welty, Hobbs, & Kulkarni 1994; Welty & Hobbs 2000). The change in 6- value affected only 
the column of Na I since the optical depth at fine center is larger for Na I than for K I. 
This resulted in a 20 % decrease in A'"(Na I) toward HD 37903 and a 10 % decrease toward 
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HD 147009; no change in A'"(K I) occurred with the change in &-value. 
Editor place figures 1-6 here. Editor place tables 1-5 here. 

2.2. Ultraviolet Data 

High-resolution ultraviolet spectra of HD 37903 and HD 200775 were available in the 
lUE NEWSIPS archive. Tabic 6 lists the short-wavelength spectral images used in our 
study. The data were extracted with the set of programs in lUERDAF. For features that 
appeared in two orders, the cleaner looking spectrum was used for further analysis. Once 
extracted, the spectra were analyzed with the IRAF package in much the same way as 
described above for the ground-based measurements. Our focus was on absorption from 
^^CO, ^^CO, C I, S I, and Ni II; sample spectra are shown in Figures 7 — 8. Since C I and 
S I represent minor constituents for their respective elements and have similar ionization 
potentials, the S I lines provided a check on the C I lines. On the other hand, Ni II is the 
dominant form of Ni in cloud envelopes; the synthesis of Ni II lines allowed for a comparison 
of 6-valucs which could elucidate the relative volumes occupied by dominant versus minor 
forms of an clement. 

Editor place table 6 here. Editor place figures 7-8 here. 

The compilation of Wx values appears in Table 7. Joseph et al. (1986) reported results 
for HD 37903 from the same lUE spectra, but with an earlier version of the archiving 
process. The overall agreement between the two sets of Wx is good; the NEWSIPS data 
used here produce spectra with higher signal to noise, as seen by the smaller uncertainties 
in Wx- We performed several checks on Wx^s measured from NEWSIPS spectra; these are 
described in the Appendix. For the features considered in our study, we are confident that 
they have been measured reliably. 

Profile synthesis was used to extract b-values and column densities from the spectra. A 
chi-squared minimization procedure was adopted for the fits. The syntheses were based on 
the following set of oscillator stengths: Morton & Noreau (1994) for ^^CO and ^^CO; Morton 
(1991) and Zsargo, Federman, & CardeUi (1997) for C I; Federman & CardeUi (1995) for 
S I; and Morton (1991) and Zsargo & Federman (1998) for Ni II. The resulting column 
densities appear in Table 8. For the synthesis of the CO bands, the excitation temperature 
was set at 4 K (Lambert et al. 1994) and the 6-value was found to be somewhat less than 1 
km s~^. For the atomic lines, the 6- value was determined by fits to several lines. There was 
no discernible difference in 6-value between neutral atoms and ions; we found b = 2.0 ± 0.5 
km s"-'^. Since CO is expected to occupy a more restricted volume than even neutral atoms. 
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a smaller 6-value for this species is appropriate. 



Editor place tables 7-8 here. 



3. Analysis and Results 



3.1. Chemical Modeling 



The chemistry of the simple molecular species, CH, C2, and CN, is relatively well 
understood (van Dishoeck & Black 1986; Federman et al. 1994) and relevant reaction rates 
are known or have been calculated. The purpose of obtaining data on several species of 
carbon-bearing molecules is to glean information about the physical conditions where these 
species reside. In particular, gas density (n), gas temperature (T), and the enhancement of 
the ultraviolet flux over the average interstellar value given by Draine (1978), can be 
extracted from the measurements. 

The analysis is based on steady-state chemical rate equations, which were successful 
in reproducing data for quiescent diffuse interstellar clouds (e.g., Federman et al. 1994) 
and the PDR toward HD 200775 (Paper I). In this analysis, the three quantities, n, T, and 
luv are inferred from observed column densities of CH, C2 and CN. However, only upper 
limits to the column densities are available for the C2 molecule toward HD 37903 and 
HD 147009 and for CN toward HD 147009, and as a consequence, the physical conditions 
in the foreground gas are not as well constrained as in other clouds (e.g.. Paper I). 

The rate equations used here arc based on our earlier work. Federman and Huntress 
(1989) developed the chemistry for C2 in diffuse clouds, while Federman, Danks, and 
Lambert (1984) and Federman and Lambert (1988) did the same for CN. The incorporation 
of NH by Federman et al. (1994) into the rate equation for CN was the result of the 
discovery of NH in cloud envelopes (Meyer & Roth 1991). The respective equations for C2 
and CN are: 



iV(C2) 



fci.x(C+)iV(CH)ntt 



and 



(1) 



G{C2) + k2x{0)n + hx{N)n 



N{CN) 



[k3xiN)N{C2) + hx(N)N{CR) + hx{C+)N(NR)a\n 
G(CN) + kex{0)n 



(2) 



In these equations, the reaction rate constants are denoted by ki (see Table 9); the 
abundance of species X relative to the proton abundance and the column density of species 
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X are given by 2;(X) and A^(X), respectively. The total proton density [n(H) + 2n(H2)] 
is given by n; the photodestruction rate of species X is G'(X); and a accounts for the 
conversion of C"*" into CO (Federman et al. 1994). The rate constants for neutral-neutral 
reactions include a factor, (T/300)°'^, to account for the velocity term in the rate constant. 
The attenuation of radiation by dust is included in the photodissociation rate, G'(X) = 
Go /u^exp(-r„^), where Go is the rate at the surface of the cloud, and is the optical 
depth at 1000 A. A wavelength of 1000 A for was chosen because it is in the middle of 
the range of wavelengths that leads to photodissociation of C2 and CN (912 A < A < 1200 
A) (Pouilly et al. 1983; Lavendy, Robbe, & Gandara 1987). The observed values of A^(CH) 
and N{C2) were used in determining the predicted values of N{C2) and A^(CN). The value 
of N(NE) was set to 0.041A^(CH) since Ar(NH) scales with A^(CH) (Meyer & Roth 1991; 
Federman et al. 1994). 

Editor place table 9 here. 

The following considerations apply to the parameters, r„^, a, and a:(X). From 
knowledge of the ratio of total to selective extinction, — A^/ E{B — V), and the 
ultraviolet extinction curve (CCM89; Buss et al. 1994), information on can be derived. 
For the direction toward HD 37903, Tuv was set at 2A„, with an somewhat larger than 
average (Buss et al. 1994). Since the extinction curve for HD 147009 is similar to that seen 
for p Oph, toward p Oph (Federman et al. 1994) was scaled by E{B — V) to obtain the 
value for the sight line to HD 147009. For optical depths greater than 2, an additional term 
of the form, a = [I + IA{tuv - 2)/5]-^ (Federman et al. 1994), was added to account for the 
trends between the predicted and observed A^(C2) and r^v caused by C+ being converted 
into CO (Federman & Huntress 1989). Fractional interstellar abundances of 1.4 x 10~^ 
(Sofia et al. 1997), 7.5 x 10"^ (Meyer, CardeUi, & Sofia 1997), and 3.2 x 10"^ (Meyer, Jura, 
& CardeUi 1998) for C+, N, and O, respectively, were adopted from measurements with 
the Hubble Space Telescope (HST). These abundances account for modest depletion of the 
elements onto interstellar grains in cloud envelopes. These values are somewhat different 
than the fractional abundances quoted in Paper 1: while the O and N abundances differ 
by 10%, the C+ abundance is almost a factor-of-2 less. The lower fractional abundance of 
carbon requires a higher density for self-consistent chemical results, mainly because C2 is 
formed through reactions involving C+ and CH. 

The column densities derived from the rate equations compare favorably with the 
observed columns for C2 and CN. The predicted column densities and inferred physical 
conditions are found in Table 10. The Table also displays predictions for the molecular 
columns toward HD 200775 based on the atomic abundances adopted here; the greatest 
effect was on the extracted density. Changing the temperature did not significantly alter 
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the predicted column densities in the chemical analysis above. The temperature used in 
this analysis was 40 K, the average value obtained from the C I excitation analysis (T ~ 
50 K) and the temperature determined from CO fractionation (T ~ 30 K). The average of 
these two temperatures is reasonable in light of the fact that C I is more extended than 
CO in the cloud. The value of luv/n is the most important parameter in PDRs (HT99); it 
controls the chemistry and structure. When n was increased in our model, had to be 
increased in order to obtain the best match to the observed column densities. The ratio 
luv/n remained unchanged (see eqn. 1 and 2). 

Editor place table 10 here. 



3.2. Atomic Analysis 

In this analysis, only the columns of Na I and K I associated with the molecular 
material (i.e., at the same velocity as the molecules CH, CN, and C2) are used. The 
atoms Na and K are predominantly singly ionized in interstellar space. Through ionization 
balance, we can analyze the neutral columns with the equation. 



iV(X) = 



Ntot{li)A{X)a{X)neac 
G{X) 



(3) 



In this equation, ^(X) is the elemental abundance, NtotiH) is the total proton column 
density, iV(H) + 2A^(H2), a(X) is the radiative recombination rate constant (Pequinot & 
Aldrovandi 1986), G'(X) is the photoionization rate for species X, and rig is the electron 
density. ^(X) includes depletion onto grains at about 25% for alkali elements (Phillips, 
Pettini & Gondhalekar 1984) and Oc incorporates the conversion of 0+ into CO in the 
molecular material (as opposed to a in the chemical analysis which is an average for the 
total amount of neutral material). This conversion is a significant effect only when the 
fraction of molecular gas is > 50% of the neutral gas. The form of ccc is a/ + (1 — /), 
where / is the II2 fraction (2A^(H2)/A^tof (H)) and a is described in §3.1. Many of the 
uncertain factors (e.g., Ue and depletion) in the above equation can be eliminated by taking 
the ratio of A^(Na I) to N{K I). This ratio can be written as 



"A^(Na I)" 
A^(K I) _ 



"A(Na) 



"Q;(Na I)" 
_ a(K I) _ 



" G{K I) " 
_G(Na I)_ 



(4) 



The appropriate value for ^(Na)/^(K) is 16 and Q;(Na l)/a{K I) is 1.05. The latter ratio 
was incorrectly given in Paper I; the change results in a 10% increase in the predicted ratio 
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of Na I and K I columns. The predicted ratio mainly depends on the photoionization rates. 
The photoionization rates for Na I and K I at the cloud surface were determined from the 
measured ionization cross sections of Hudson & Carter (1967a) and of Hudson & Carter 
(1965, 1967b), Marr & Creek (1968) and Sandner et al. (1981), respectively. 

To complete the atomic analysis, details about the extinction curve are needed. The 
main difference for the two sight lines studied here involves the effect of extinction from 
interstellar grains. The pair method (WBS84) was used to derive the extinction curve for 
the central star in the nebula vdB 102, HD 147009. Extinction curves for two other stars, 
HD 37023 and HD 37903, were computed for comparison with extinction curves derived 
by Fitzpatrick & Massa (1990 - hereafter FM90). Stellar spectra of HD 31726, HD 37023, 
HD 37903, HD 55857, HD 64802, HD 147009 and C Aql were obtained from the lUE 
Archive; the filenames are hsted in Table 6. The unreddened standard stars used in FM90 
are used in this work, with the exception of C Aql (Witt, Bohhn, & Stecher 1981), in order 
to obtain the most consistent results. The extinction curves for HD 37023 and HD 37903 
derived by FM90 are similar to those derived by us, demonstrating that the extinction 
curve for HD 147009 derived here is reliable. 

The extinction curves for HD 37903 and HD 147009 were found to be anomalous when 
compared to the one representing typical Galactic extinction, that for Per (Savage & 
Mathis 1979). For this analysis, the extinction curve for C Per was scaled to the same 
reddening parameter of the observed star. Figure 9 displays the resulting extinction curves. 
Since the photoionization rates were determined for ^ Per (Paper 1), the difference between 
extinction toward the reddened line of sight compared to the average is used to derive an 
attenuated photoionization rate for the reddened star. Changes in extinction introduce 
a multiplicative factor to the photoionization rates of the form 10^. For the wavelengths 
where the curve for C Per shows more (less) extinction, the exponent /? is positive (negative). 
Since the two extinction curves can intersect, for some wavelength intervals the extinction 
can be enhanced and diminished for others. The total rate is the sum of the rates for each 
wavelength interval. 

Editor place figure 9 here. 

3.3. Atomic Results 

The atomic spectra for Ca II H and K and Na I D exhibit 6 velocity components along 
the line of sight toward HD 37903 (see Figure 1). There is no molecular absorption seen 
in 5 of them. The atomic gas not associated with molecular absorption is attributed to 
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Orion's Cloak. Orion's Cloak is a rapidly (35 < v < 120 km s~^) expanding shell of gas 
centered on the Orion OBI association (Cowie, Songaila & York 1979). While the 5 velocity 
components toward HD 37903 have lower velocities (< 30 km s~^) than is seen in Orion's 
Cloak, this velocity structure is observed toward many other stars in Orion (Cowie et al. 
1979; Welty ct al. 1994; Welty et al. 1996) showing the high velocity gas of Orion's Cloak. 
Our atomic analysis involves only absorption associated with the molecular material in the 
nebula. 

The atomic observations yield information on the photoionization rates present in the 
gas (see eqn. 4) because eqn. 4 is not dependent on the physical conditions. Toward ( Per, 
an observed ratio Na I/K I of 100 ± 6 was obtained. For the predicted ratio a refined 
analysis was performed. The difference in grain attenuation between 2400 and 2900 A {AAx 
— 0.2) was taken into account. This resulted in a predicted ratio of 67. Toward HD 200775, 
the observed ratio is 54 ± 9, which is in excellent agreement with the predicted ratio of 
54 when the effects of the extinction curve are applied. (Note that the sum of all velocity 
components of Na I and K I are used for this sight fine, see §4.4.) 

The observed column densities were (1.1 ± 0.2) x 10^^ cm~^ and 1.1 x 10^^ cm~^ for 
Na I and K I toward HD 37903, giving an observed ratio of 100 ± 18. The inferred Na I and 
K I photoionization rates are 1.2 x 10^^^ s^^ and 4.4 x 10~^^ s~^, which yield a predicted 
ratio of 62. For the line of sight to HD 147009, the observed column densities of (1.5 ± 0.3) 
X 10^^ cm~^ and (2.3 ± 0.1) x 10^^ cm^^ for Na I and K I, respectively, were determined. 
These columns produce an observed ratio of 65 ± 13. The Na I and K I photoionization 
rates for this gas are 2.9 x 10~^° s~^ and 9.7 x 10^^° s~^, respectively. These values lead 
to a predicted ratio of 56. As expected the correspondence between observed and predicted 
N{N& l)/N{K I) ratios is good. 

Editor place table 11 here. 

Another useful way to study the observed atomic abundances is to compare, for a 
specific atom, the column density obtained toward the program stars to the column density 
obtained toward ( Per. In this way we check for potential differences in depletion onto 
grains and He (see eqn. 3). Curves of growth with b = 1.5 km s~^ for the weak Na I 
doublet at 3302 A (Crutcher 1975) as well as the strong D-lines at AA 5889, 5895 (Hobbs 
1974a, 1974b; Welty et al. 1994; Andersson, private communication) and the weak K I 
doublet at AA 4044, 4047 (Knauth et al. 2000) yielded weighted average column densities 
for Na I of (7.3 ± 0.3) x lO^^ ^ni-^ and (7.3 ± 0.3) x 10^^ cm-^ for K I toward C Per. 
Our precise column density measurements of K I AA4044, 4047 yield similar results to the 
column density obtained from the K I A7699 line. The photoionization rates of Na I and 
K I toward ( Per are 4.5 x 10~^^ s~^ and 1.8 x 10~^^ s~^ (Paper I). For this analysis we 



-14- 



require the total proton density toward each star. The proton density is (1.6 ± 0.3) x 10^^ 
cm"^ toward ( Per (Bohhn, Savage, & Drake 1978) and is (3.6 ± 0.7) x 10^^ cm~^ toward 
HD 37903 (Buss et al. 1994). Toward HD 147009, the total proton column density was 
determined from several sources. From radio observations of CH, the column density of H2 
was found to be 2.5 x 10^° cm~^ (Vallcc 1987) using the relationship derived by Federman 
(1982) and Banks, Federman, & Lambert (1984). A value of N{}i I) = 1.1 x lO^^ cm'^ 
was determined from the 21-cm results of Cappa de Nicolau & Poppel (1986), indicating 
Ntot(ii) ~ 1-6 X 10^^ cm~^. Also from 21-cm observations, de Geus & Burton (1991) derived 
Ar(H I) = 1.8 X 10^^ cm-^ which yields Ntot{^) = 2.3 x 10^^ cm-^. The average proton 
column density toward HD 147009 is therefore (1.9 ± 0.4) x 10^^ cm~^. In our analysis, we 
adopted a 20 % error in Ntoti^) for all stars except HD 200775 where the observed precision 
was higher. 

Most of the total proton density is assumed to lie in front of HD 147009, since most 
of A'"(H I) resides in one component (Cappa de Nicolau & Poppel 1986) and the star is 
at high galactic latitude, h° = 20.894 (GZL89). The relationship of Bohlin et al. (1978), 
Ntot(ii)/ E{B — V) = 5.8 X 10^^ atoms cm~^ mag~^, supports this assumption. This 
relationship yields Ntot{ii) = 2.2 x 10^^ cm^^, which is similar to the column density 
determined above. Therefore it is reasonable to assume all of the neutral hydrogen resides 
in front of the star. 

Estimates for gas density, temperature and UV field strength are needed to complete 
this analysis. The values of n, T, and luv toward ^ Per (Federman et al. 1994) were 
re-derived with the new values for the fractional abundances of C"*", N, and O, resulting in 
1200 cm-^ 40 K, and 1. For the gas toward HD 37903 and HD 147009, the adopted values 
of n and T are from the chemical analysis above for luv — 1 (see Table 10). Here we set 
rie — n a:(C"^), since the electrons mostly come from ionized carbon. The predicted ratios 
toward HD 200775 are 0.05 ± 0.01 and 0.06 ± 0.01 for Na I and K I, respectively. The 
observed ratios are 0.08 ± 0.01 and 0.15 ± 0.02 for Na I and K I, respectively. (These ratios 
are based on our refined analysis and differ somewhat from the ratios in Paper 1.) For the 
case of HD 37903, the effect of needed to be included since the fraction of molecular gas 
is > 50 %. This resulted in a predicted ratio A^(Na 1)hd 37903/MNa l)^ Per of 0.23 ± 0.06 
and 0.21 ± 0.06 for K I. The observed ratios are 0.15 ± 0.03 for Na I and 0.15 ± 0.01 for 
K I. As for the PDR toward HD 147009, the analysis gives a ratio of 0.01 ± 0.01 for both 
ratios. The observed column density ratios are 0.21 ± 0.04 for Na I and 0.32 ± 0.02 for K I. 

The results of our analysis are summarized in Table 11. It can be seen from the ratios 
of Na I and especially K I that a modest change in rig and/or depletion onto grains is 
needed to bring the predictions and observations for the foreground PDR of NGC 7023 
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into agreement. This change from Paper I is caused by the decrease in a;(C+). Since the 
observed and predicted results for Na I and K I toward HD 37903 with respect to ( Per 
agree within the errors, our estimates for and depletion for the foreground material 
toward HD 37903 are reasonable. The situation is more extreme for the foreground PDR 
of vdB 102. A factor-of-20 increase in Ue is needed to bring the predictions into agreement 
with the observations. While the fiat FUV extinction curve for HD 147009 allows more 
radiation to permeate the nebula, rig is limited by the elemental abundance of carbon. 
Because HD 147009 is an AO V star, enhanced ionization in an H II region is not expected, 
unless the enhanced ionization arises from the combined effects of other, earlier type, stars 
in Sco 0B2. 



3.4. UV Results 

The lU E data on C I and CO provide additional constraints on the physical conditions 
for the foreground PDR. The relative populations of fine structure levels in the electronic 
ground state of C I arises from a combination of coUisional excitation and de-excitation, 
radiative decay, and in enhanced radiation fields, UV pumping. Analysis of the observed 

distribution via statistical equilibrium yields pressure {ritot T) and sometimes J„^. Here, ritot 
is the density of collision partners, n(H) + n(H2). If the kinetic temperature is known from 
another diagnostic, ritot can be obtained and compared with n [n(H) + 2 ?^(H2)] inferred 
from the chemistry. The ratio of A^(-'^^CO)/A^("'^'^CO) is affected by chemical fractionation. 
When the gas is cold (T < 40 K) and is permeated by a strong UV radiation field so that 
the abundance of C+ is high, the ratio of columns is less than the ratio of ambient ^^C/^^C 
(about 60) because ^^CO has the lower zero-point energy (cf. Watson, Anicich, & Huntress 
1976). If the gas in not especially cold, but still affected by UV photons, selective isotopic 
photodissociation enhances A'"(^^CO) relative to iV(^^CO) (e.g.. Bally & Langcr 1982; van 
Dishoeck & Black 1988). The enhancement is the result of ^^CO self-shielding itself from 
photodissociation because its lines are more optically thick. 

The analysis of C I excitation described by Lambert et al. (1994) forms the basis of 
our results. In particular, their collisional rate constants and radiative rates are adopted 
here. The fractional abundances of collisional partners arc based on the observations of 
Buss ct al. (1994) - sec Table 12, with the added assumption that 10% of the gas is He. 
The analysis also assumes that absorption from each fine structure level, J, arises from the 
same parcel of gas. Figure 10 shows the results of our calculations. Here we plot the density 
and temperature consistent with the observed column density ratios, A^(J = 1)/A^(J = 0) 
— dashed curves — and A^(J = 2)/A^(J = 0) — dash-dotted curves. The upper and lower 
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bounds are based on ± 1-a errors in column density. For HD 37903, results for of 3 are 

indistinguishable from those with 1^^ of 1. The strongest constraint comes from N{2)/N(0): 
ritot lies between 300 and 500 cm~^ and T about 50 K. Based on the observed columns of 
H and H2, the inferred value for n is 400 to 700 cm^^. As for the foreground PDR toward 
HD 200775, the lower limits from A^(2)/A^(0) and the upper limits from A^(1)/A^(0) indicate 
that ritot is about 700 cm~^ (or n is ~ 850 cm~^) for T between 40 and 60 K. Analysis of 
C2 excitation from the data in Paper I is consistent with the C I results if /jj. is 3. For 
C2, optical pumping to an excited electronic state involves absorption of infrared radiation; 
hence the use of lir- [Again, the curves for I^y of 3 are essentially those for I^v of 1, except 
for the lower range allowed by A'"(2)/A'"(0).] 

Editor place table 12 here. Editor place figure 10 here. 

The inferred values for n agree with the chemical results for the two foreground PDRs. 
The gas containing C I toward HD 37903 appears to have a density most consistent with 
the chemical results for 1^^ — 1. Since C I may occupy a larger volume of material than 
the molecules C2 and CN, the inferred average density from C I excitation could be smaller 
by about a factor of 2 compared to the density in CN-rich gas. This seems to be the case 
for C Oph (Lambert et al. 1994). In other words, the chemical results for of 1 — 3 
are still consistent with the results from C I. The situation is different for the gas toward 
HD 200775. Here, the C I results are most consistent with the chemical results for I^v of 3. 
Since this comparison also applies to the results for C2 with lir of 3, the strengths of the 
UV and IR fields seem to scale together for the foreground PDR. 

The observations of HD 200775 reveal that N{^^CO)/N{^^CO) is about 20. An 
enhancement in A^(^^CO) arises when the gas is cold, yet contains an appreciable fraction 
of C^. Lambert et al. (1994) described the amount of fractionation in terms of the quantity 
Fi3, the column density ratio divided by the ambient ^^C/^^C ratio. When the ratio of 
columns is smaller than the ambient ratio, F13 equals exp(-35/r). Since the gas toward 
HD 200775 has F13 of ^ 0.33, the CO observations indicate a temperature of 30 — 35 K. 
The need for C+ suggests that the UV radiation field is somewhat higher than in a typical 
diffuse cloud, consistent with the results described above. Since ^^CO was not detected 
toward HD 37903, no meaningful analysis is possible. 



4. Discussion 
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4.1. Structure of PDRs 

Since /„^,/n determines the structure of PDR's (HT99), we can estimate the fraction of 
foreground material in the PDR. We use eqn. 7 from Hollenbach & Tielens (1999), which 
approximately relates the formation rate of H2 to the photodissociation rate, for Ntot{ii) 
< 10^^ cm~^. Due to its approximate nature, we assume that the use of this equation for 
lines of sight with Ntot(ii) of a few times 10^^ cm~^ to be valid as well. This steady state 
expression is used to reproduce the column of H2, assuming the observed H2 is completely 
within the PDR: 

In this equation. No is the column density at which self-shielding becomes important {Ng 
= 1 X 10^^ cm~^), is the formation rate of H2 on grains {'Jhz = 3 x 10"^''' cm^ s~^), 
and /diss(O) is the photodissociation rate (4 x 10~^^/„t, s~^). Table 12 shows that the 
predicted values for luv are 30 to 60 times greater than the values derived from the chemical 
analysis. The luv determined from the chemistry must be the same as that determined from 
the above equation for H2; these results suggest that only a small fraction (~ 3%) of the 
foreground material is actually associated with the PDR in front of the star. 

The linear extent of the foreground PDR (Ltr) can be compared to that in the 
background molecular cloud. We adopted the expression for Ltr of Federman, Glassgold, & 
Kwan (1979): 



2.93 x 10^^n-^-\5.6 x IQ-^HuvY'^ 
= cm. (6) 



For the foreground PDR of NGC 2023, this equation yields an extent of 0.002 — 0.007 pc 
for densities of 500 — 1500 cm~^ and respective values for of 1 — 3. For the background 
cloud with n of 10^ cm~'^ and luv of 10"^ (Rouan et al. 1997), the extent is 0.0086 pc. For 
NGC 7023, the range in extent is 0.008 — 0.03 pc when n of 300 — 900 cm~^ and respective 
luv of 1 — 3 are used for the foreground cloud. Utihzing the density (10*^ cm~^) and of 
10^ from Rogers, Heyer & Dewdney (1995) for the background cloud indicates a length scale 
of 0.086 pc. The extent of the PDR into the molecular cloud behind the stars is comparable 
to its extent in the lower density foreground gas. For vdB 102 information is only available 
for the foreground gas; an extent of 0.004 — 0.01 pc is found from densities and luv of 400 
— 1200 cm~^ and 1 — 3, respectively. 



Prom the column of neutral hydrogen and the density of the cloud, one can estimate 
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the size of the neutral cloud in front of the star, assuming all the hydrogen is associated 
with the cloud. Toward HD 37903, a neutral cloud size of 0.8 — 2.3 pc is derived. A cloud 
size of 0.5 - 1.5 pc toward HD 147009 and toward HD 200775 a cloud size of 0.8 - 2.5 pc 
are obtained. These simple calculations show that the PDR does indeed reside in a small 
fraction (~ 1 %) of the neutral cloud. 

4.2. NGC 2023 

In order to arrive at a consistent model of PDR's associated with these reflection 
nebulae, it is necessary to compare the results obtained here to the physical conditions 
inferred from diagnostics at other wavelengths of the background molecular cloud. As noted 
in §1.1, there have been many studies of NGC 2023 whose focus was the molecular cloud 
behind HD 37903. Here, we highlight those studies that relate to our measurements. 

Infrared H2 observations (Martini et al. 1999 and references therein) of this nebula 
show that there is a neutral shell of molecular material 0.13 — 0.20 pc from the star. 
However, this neutral shell is not seen in the H I data. The H I emission is observed to be 
centered on HD 37903 (Lebron & Rodriguez 1997). Clumps of material have been detected 
in the nebula on less than parsec size scales (Martini et al. 1999). This implies that there 
are two components within the nebula. Steiman-Cameron et al. (1997) determined that the 
densities for the low and high density components are 750 cm~^ and 10^ cm""^, respectively, 
with a temperature range of 250 K < T < 750 K and a radiation field of 1.5 x 10^ times 
the average interstellar radiation field. For the gas directly behind HD 37903, Wyrowski 
et al. (1997) found a density of 3 x 10*^ cm~^, a radiation field ~ 500 times the average 
interstellar radiation field (ISRF) and a temperature of 125 K. 

Martini et al. (1999) measured the near-infrared emission from H2 in two regions, 60" 
south and 160" north of HD 37903. Their results show that these regions are of relatively 
high density, n ~ 10^ — 10^ cm^'', with a ratio of radiation field to density /^/n ~ 0.1 — 
0.01 cm^. These measurements are in agreement with other studies which find high density 
clumps (10^ — 10^ cm~^) in a slightly less dense interclump medium (~ 10*^ — 10^ cm^'^) 
with a enhancement of the UV radiation field over the average ISRF of 10^ — 10^ for 
NGC 2023 (Martini et al. 1999, Takami et al. 2000). 

Our chemical analysis reveals that the radiation field is weaker in the foreground 
material, ~ 1 — 3, than the values determined in previous work for the molecular cloud 
behind the star (Field et al. 1994, 1998; Fuente et al. 1995; Steiman-Cameron et al. 1997; 
Wyrowski et al. 1997). The density of the foreground gas has a range of a few times 10^ 



-19- 



to 10^ cm~^, a factor of 10 to 1000 less dense than the gas in the molecular cloud behind 
HD 37903. The temperature of the foreground gas is a factor-of-2 less than the material 
behind the star. These observations of the foreground material indicating lower densities 
support the idea that HD 37903 formed near the edge of the molecular cloud and rapidly 
dispersed the material at the edge of the cloud. NGC 2023 could be a face-on example of 
star formation by the blister model (Zuckerman 1973; Israel 1978). 

Burton ct al. (1998) indicate clump densities of ~ lO"' cm^'^ and a far ultraviolet 
radiation field about 10'^ times the ambient interstellar value at a projected distance of 0.18 
pc. Buss ct al. (1994) measured a stellar flux at 1050 A of 1.5 x 10"^^ erg cm~^ s^^ A^^. 
Using an average interstellar flux of 2.65 x 10^^ ergs cm^^ s^^ A.~^, = 2.9, and a 
distance to the star of 470 pc (Ferryman et al. 1997), we find a distance to the foreground 
gas of 2.5 — 4.5 pc for our values of between 1 and 3. Therefore the radiation field at 
a distance of 0.1 — 0.2 pc from the star would be 500 — 2000 times the average radiation 
field, which is consistent with Burton et al.'s (1998) value of ~ 1000 (usually denoted 
by Go) at a distance 0.18 pc from the star. 

Prom maps of CN J = 1 ^ (113.5 GHz) and HON J = 1 ^ (88.6 GHz) emission, 
Fuente et al. (1995) detected an enhanced amount of CN in the molecular shell in 
NGC 2023. The CN/HCN ratio inside the H2 shell was found to be 25 times larger than 
outside. Outside the shell this ratio is close to the interstellar value for molecular clouds (~ 
2). HON photodissociation plays an important role because it produces CN and because the 
photodissociation rate of HCN is greater than that of CN. While we previously found an 
enhancement in the abundance of CN toward HD 200775 (Paper I), consistent with radio 
observations of NGC 7023 (Fuente et al. 1993), our chemical results for the foreground 
FDR of NGC 2023 show no such enhancement. Most of the CN is in the envelope of the 
molecular cloud, not in the foreground FDR. 



4.3. vdB 102 

The chemical analysis for HD 147009 matches well with the observed upper limits 
for C2 and CN, but the physical conditions derived in our chemical analysis could not 
be constrained well. Unfortunately, little is known about the physical conditions in this 
nebula. Kutner et al. (1980) mapped CO emission in refiection nebulae and determined 
that the peak CO emission is offset, about 3.3 arcmin to the north-east, from the star. 
From radio observations of CH in the cloud, a column density 7V(CH) = (1.0 ± 0.3) x 10^^ 
cm~^ was obtained (Vallee 1987). This is comparable to the column density obtained from 
our measurements, A'^(CH) = (1.5 ± 0.1) x 10^^ cm~^. The peak CH emission reported 



-20- 



by Vallee is at the same offset as Kutner's CO observations. Since there is a proportional 
relationship between iV(CH) and A^(H2) (Federman 1982; Banks et al. 1984; Mattila 1986), 
Vallee inferred that A'(H2) = 2.5 x 10^° cm~^. Vallee used a length of 0.13 pc to convert 
column density to space density. This length was determined from the fact that the CH 
abundance extends over 0.3 magnitudes in blue absorption, As- There could easily be a 
factor-of-2 error in this estimate. With this length, Vallee obtained n(H2) = 600 cm~^ and 
then converted 71(^2) into a total proton density oi n — 1300 cm~^ with the relation n — 
2.2 n(H2), assuming 10% of the gas is He. Vallee's density, which is similar to ours for 
= 3, may not represent the molecular cloud behind the star as he reported, since his CH 
abundance is comparable to ours. The molecular material may in fact reside in front of 
HD 147009. Moreover, this density is similar to densities found in the diffuse component in 
the reflection nebulae NGC 2023 (this work) and NGC 7023 (Paper I and this work). 

4.4. Optical Results for Three Nebulae 

The reflection nebulae NGC 2023 and NGC 7023 are relatively bright extended objects, 
located at the edges of their respective molecular clouds. For these reasons they have been 
extensively studied at a variety of wavelengths. The two nebulae have been shown to have 
clumpy, filamentary structure from infrared observations. The central stars, HD 37903 and 
HD 200775, of these nebulae have similar spectral types, B1.5 V and B3e V, and similar 
distances of 470 pc and 430 pc (Ferryman ct al. 1997), respectively. The peak infrared 
emission is offset from the stars in both nebulae. As for vdB 102, the illuminating source, 
HD 147009, is of later spectral type (AO V) and is closer to the Sun at a distance of 160 pc 
(Ferryman et al. 1997). The extent of filamentary structure is unknown. 

The physical conditions obtained from optical spectra toward the illuminating stars 
in the reflection nebulae NGC 2023 and vdB 102 are similar to those obtained from an 
earher study of NGC 7023 (Faper I) and updated here. The physical conditions n, T, and 
luv inferred for NGC 7023 are 300 — 900 cm~^, 40 K, and 1 — 3, respectively, once the 
effects of a lower a;(C+) are included (Table 10). These physical conditions are similar to 
those obtained from grain scattering observations of NGC 7023 (Walker ct al. 1980; Witt & 
Cottrell 1980a, 1980b; Witt et al. 1982). The physical conditions found toward NGC 2023 
are 500 - 1500 cm^^ 40 K, and 1 - 3 and toward vdB 102 400 - 1200 cm'^, 40 K, and 
1 — 3, for n, T, and respectively. Therefore, the material in front of the stars in the 
reflection nebulae NGC 2023, vdB 102, and NGC 7023 is relatively diffuse. 

In Faper 1, an enhancement in the CN abundance over that described in eqn. 2 
was (and is still needed) to explain the observed CN abundance in the foreground FDR 
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of NGC 7023. This enhancement of CN in the foreground material of NGC 7023 is the 
result of a contribution from dark cloud chemistry (Paper I) — HCN — ^ CN through 
photodissociation. If we really are seeing an effect of dark cloud chemistry from the 
molecular cloud behind the star, the foreground PDR may be closer to the star than our 
estimates suggest, even for an enhancement by a factor-of-3 in the radiation field. This 
enhancement is not detected toward either HD 37903 or HD 147009, unless the amount of 
C2 in the gas is much less than our current upper limit. This could be a consequence of the 
different distances between each star and its foreground gas. 

In the case of NGC 7023, an hourglass shaped cavity has been detected (Lcmairc et al. 
1996; Fuente et al. 1998) around the star. CO radio maps reveal (Gerin et al. 1998) that 
the ^^CO emission profiles are relatively smooth and centered at 2.5 km s~^, while the ^^CO 
profiles indicate more complicated structure. This structure shows a gradual shift in the 
velocity of the gas from 2.5 km s~^ to 1.5 km s~^, which is the velocity of the foreground 
molecular material reported in Paper I. This suggests that there should be a diffuse ^^CO 
component in the cavity or the foreground material. We associate this difi^use component 
with the CO seen in UV absorption. Lemairc et al. (1999) report redshiftcd H I with 
respect to II2 in the molecular cloud behind the star. This trend is also detected in atomic 
absorption (Paper I) which is blueshifted compared to the foreground molecular material. 
This is consistent with an expanding PDR around the star (Gerin et al. 1998; Fuente et al. 
1998). 

The analysis of atomic data provides insight into the processes taking place in the 
foreground gas. Using the extinction curve for the line of sight in our model, the observed 
ratio of Na I to K I toward NGC 7023 and vdB 102 matched the predicted ratio. In the 
cases for ( Per and NGC 2023, the observed ratio was about twice the predicted value. 
Adjustments to the photoionization rates would bring the predicted ratios into better 
agreement with observations. In particular, a finer grid for computing the differences in 
extinction relative to the extinction for C Per may be more appropriate; we only considered 
average differences over wavelength regions where the curves intersected. In the comparison 
of the column density of each species toward the reflection nebulae and ( Per, a small 
change in rig or an extra source of depletion is needed toward NGC 7023 to match the 
observed and predicted results. Toward HD 147009, a fairly substantial modification is 
required. A significant increase in Ug, possibly provided by protons in diffuse ionized gas 
created by earlier type stars in Sco 0B2, is needed to account for the differences found 
between the nebula vdB 102 and gas toward ( Per. 

The analysis of ultraviolet data provides additional pieces to the picture for the 
gas associated with NGC 2023 and NGC 7023. It appears that the foreground PDR in 
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NGC 7023 is more homogeneous (because C I excitation, C2 excitation, CO fractionation, 
and chemical analysis yield very similar physical conditions) and is permeated by a slightly 
stronger radiation field compared to that for NGC 2023. 

5. Conclusions 

We presented the first chemical results for foreground material associated with the 
PDRs of NGC 2023 and vdB 102. The foreground material in NGC 2023 is more diffuse 
than the molecular cloud behind HD 37903 and is subjected to a less intense radiation 
field. The results are similar to what we found for NGC 7023 (Paper I). As for vdB 102, 
our analyses are the most comprehensive to date and will be the basis for other studies of 
nebulae illuminated by A-type stars. 

Future studies should incorporate details that extend the treatment presented here. 
The effects of a clumpy distribution of material and a more extensive set of chemical 
reactions in the framework of a time-dependent model need to be considered. These models 
have to include the "anomalous" extinction law found for both stars in this work (sec also 
FM90; WBS84). The consequences of NGC 2023 and vdB 102 being located in active 
star-forming regions is another area for future analysis. Other absorption measurements of 
species detected at ultraviolet wavelengths, such as H2, are needed to place more stringent 
constraints on physical conditions. Observations with the Far Ultraviolet Spectroscopic 
Explorer will soon be obtained. This combination of visible and UV data, on the one hand, 
and more sophisticated modeling efforts, on the other, will provide a clearer description of 
the interaction between a star's radiation field and the surrounding molecular cloud. 
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The lUE archive is now comprised of NEWSIPS data, which yields spectra with 
superior signal to noise. Techniques were introduced to perform geometric and photometric 
corrections, to improve the ripple correction, and to determine the image background (e.g., 
Nichols 1998). We examined the reliability of our data in two ways. First, we compared our 
results for another star, 20 Aql, with those of Hanson, Snow, & Black (1992) based on an 
earlier form of the archive. Then we compared Wx^s of C I lines seen in the lUE spectra 
of X Per with measurements from high-resolution spectra acquired with ECH-A of the 
Goddard High Resolution Spectrograph on the Hubble Space Telescope. Our comparisons 
show that NEWSIPS data provide a consistent, rehable set of Wx^s with improved signal to 
noise for interstellar studies, thereby confirming the analyses performed by Gonzalez-Riestra 
et al. (2000). 

Tables 13 to 15 provide the basis for comparison. Table 13 shows our results for 
interstellar gas toward 20 Aql based on NEWSIPS data and the results of Hanson et al. 
(1992). The agreement in Wx for values from less than 10 to 150 mA is excellent, without 
any discrepancies beyond the mutual uncertainties. The uncertainties associated with the 
NEWSIPS spectra also tend to be smaller. Sample spectra appear in Figure 11, and Table 
14 displays the column densities derived from the fits. In Table 15, NEWSIPS results 
for C I toward X Per are compared with measurements from ECH-A spectra. Again, the 
correspondence is excellent. Syntheses of the ECH-A data give a measure of the precision 
in column density attainable from NEWSIPS data. Since ^-values can be determined from 
NEWSIPS data to about 0.5 km s~^ and 0.1 km from ECH-A data, column densities 
from relatively strong hues in NEWSIPS data are accurate to about 50%. This confirms 
the measure of uncertainty given to atomic column densities in Table 8. 

For completeness. Tables 16 to 18 indicate a synthesis of results from ground-based 
observations of 20 Aql, incorporating data acquired during observing runs described here. 
This compilation includes the highest-resolution spectra taken of 20 Aql. Our spectra are 
displayed in Figures 12-14. Analyses, similar to those for NGC 2023 and vdB 102, were 
applied to the data. The density, temperature and radiation field arc 1400 cm^^, 50 K and 
1; the factor-of-2 increase from Fcdcrman, Strom, & Good (1991) in the derived density is 
due to the lower fractional abundance of carbon. With the exception of the density, the 
results, derived from the chemical analysis, are much like the results of Federman et al. 
(1991) and Hanson et al. (1992). 

The atomic analysis performed above is repeated for data toward 20 Aquilae. The 
inferred Na 1 and K I photoionization rates are 3.4 x 10^^'^ s^^ and 1.1 x 10^^^ s^^, which 
yield a predicted ratio of 54. The observed column densities of (1.3 ±0.1) x 10^^ cm~^ 
and (1.4 ± 0.1) x 10^^ cm~^ for Na I and K I, respectively, produce an observed ratio of 93 
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± 10. The correspondence between observed and predicted iV(Na l)/A'"(K I) ratio is about 
two to one, similar to that found for gas toward ( Per and HD 37903. 

For the comparison of the columns Na I and K I of 20 Aquilae to those of ( Per, the 
total proton density toward 20 Aquilae is needed. The proton density is (3.0 ± 0.6) x 
10^^ cm~^ (Hanson et al. 1992) with an assumed 20 % error (see text above). This results 

in a predicted ratio A^(Na 1)20 AquUae 

/N{Na I)^ Per of 29 ± 8 and 36 ± 10 for K I. The 
observed ratios are 18 ± 2 for Na I and 19 ± 2 for K I. The effect of was not needed in 
this analysis (see text above). Formally, the predicted and observed ratios are consistent 
with (at most) a modest change in Ue and/or depletion to bring the predictions into better 
agreement with the observations, much like the results for the FDR toward HD 200775. 

Editor place figures 11-14 here. Editor place tables 13 - 18 here. 
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TABLE 1 



Compilation of Ground-based Results for HD 37903 







McD (Moderate Res)" 


McD (High Res)" 


Other 




Species 


Wavelength 




^LSR 




fLSR 


Wx 


"LSR 




(A) 


(mA) 


(km s-i) 


(mA) 


(km s-i) 


(mA) 


(km s-i) 


CH 


4snn SIS 


7.0 ± 0.1 


11.7 


5.9 ± 0.4 


10.3 


< 7.8'' 
















< 12.5'^ 
















6.4 ± 0.3'* 


10.0 




3878.774 


< 0.6 
















1.1 ± 0.3 


11.8 






< 14.0' 






SSQfl 217 














CH+ 


4232.548 


10.7 ± 0.1 


11.8 


9.8 ± 0.3 


9.7 


10.7 ± 2.4^^ 


7.4 














10 8 -1- 5"* 


9.1 




3Q'i7 fiQ2 


5.2 ± 0.2 


10.7 






< 16.0'^ 
















5.3 ± 0.8'' 


10.1 


CN 




0.7 ± 0.3 


9.3 


< 2.5 










3874.61 


2.2 ± 0.4 


10.1 


3.7 ± 0.8 


9.1 


< 12.0'' 






3875 7fi 

(JO 1 U ■ 1 <J 


< 1.3 




< 2.5 








Ca II 




3.2 ± 0.3 


-22.9 














8.8 ± 0.4 


-13.0 














1 n Q -I- n c: 
lU.o it U.o 


-l.D 














18.0 ± 0.3 


6.2 






19.7 ± 5.0'* 








9.3 ± 0.3 


11.5 














O.O ± U.Z 


1 n 










Ca II 


3933.663 


12.8 ± 0.6 


-23.0 


11.8 ± 0.9 


-25.4 










20.7 ± 0.6 


-14.0 


22.8 ± 1.5 


-15.5 










20.3 ± 0.6 


-3.3 


17.8 ± 1.5 


-4.8 










33.6 ± 0.5 


5.9 


31.3 ± 0.7 


4.3 










15.4 ± 0.4 


11.2 


16.3 ± 1.1 


9.7 










7.5 ± 0.5 


18.8 


6.7 ± 1.3 


17.3 






Cal 


4226.728 


1.3 ± 0.4 


8.8 


< 11 




< 10.0'' 
















< 1.5'' 




Na I 


5895.924 


2.1 ± 3.0 


-20.7 














18.7 ± 3.0 


-10.5 














17.6 ± 3.0 


-4.3 














21.8 ± 3.3 


3.9 














140.6 ± 3.3 


10.0 














4.7 ± 2.2 


19.7 










K I 


7698.974 










18.8' 


10 



*2dcoude. 

'^Federman et al. 1994. 

''Allen et al. 1994. 

•^Gredel et al. 1993. 

"Welty, private communication. 



TABLE 2 



Compilation of Ground-based Results for HD 147009 







McD (Moderate Res)" 


McD (High Res)" 


McD (High Res)' 


CTIO (Moderate Res) 


Other 


Species 


Wavelength 

(A) 


Wx 
(mA) 


^Km S j 


Wx 
(mA) 


^Km s j 


Wx Vl^s-R 

(mA) (km s-i) 


Wx 
(mA) 


^Km S j 


Wx 
(mA) 


CH 


4300.313 
3878.774 
3886.409 
3890.217 


10.4 ± 0.5 

< 4.5 

< 7.5 

< 10.8 


4.3 


11.6 ± 0.8 


4.0 




12.3 ± 1.5 


-0.9 


8.8 ± 2.8^^ 


CH+ 


4232.548 
3957.692 


18.0 ± 0.7 
8.5 ± 0.9 


5.1 
2.4 


19.6 ± 1.4 


1.6 




19.0 ± 1.4 


1.8 




CN 


3874.00 
3874.61 
3875.76 


< 3.4 




< 2.4 






< 4.4 






Ca II 


3968.468 
3933.663 


13.6 ± 2.9 
29.2 ± 1.6 


1.0 
2.1 


22.4 ± 2.3 


1.4 




20.1 ± 3.2 


-0.9 




Cal 


4226.728 


< 2.6 




< 0.9 






< 2.2 






Na I 


5895.924 


111.3 ± 5.4 


2.5 






137.0 ± 4.0 1.2 








K I 


7698.974 


42.3 ± 4.0 


2.7 






33.0 ± 1.0 1.6 









*2dcoude. 

''6 foot camera. 

■^Federman et al. 1994. 



-31 - 



TABLE 3 
C2 RESULTS 

HP 37903 HP 147009 

Line Wavelength Wx N" Wx N" 





(A) 


(mA) 




(cm ^) 


(mA) 




(cm-2) 


R(0) 


8757.69 


< 


1.5 


< 


2.2 


X 10^2 


< 


2.9 


< 


4.4 X 10^2 


P(2) 


8766.03 


< 


1.2 


< 


1.8 


X 1013 


< 


2.1 


< 


3.1 X 10^3 


Q(2) 


8761.19 


< 


0.6 


< 


1.8 


X IO12 


< 


2.2 


< 


6.6 X 10^^ 


R(2) 


8753.95 


< 


0.6 


< 


2.3 


X IO12 


< 


3.2 


< 


1.2 X 10" 


P(4) 


8773.43 


< 


0.6 


< 


5.6 


X 10^^ 


< 


2.1 


< 


1.9 X 10" 


Q(4) 


8763.75 


< 


0.8 


< 


2.4 


X IO12 


< 


2.1 


< 


6.3 X 10^^ 


R(4) 


8751.68 


< 


0.4 


< 


2.0 


X IO12 


< 


1.7 


< 


7.8 X 10^^ 


P(6) 


8782.31 


< 


1.2 


< 


9.2 


X W 


< 


3.6 


< 


2.8 X 10" 


Q(6) 


8767.76 


< 


1.0 


< 


3.0 


X IO12 


< 


2.7 


< 


8.1 X 10^^ 


R(6) 


8750.85 


< 


0.4 


< 


2.1 


X IO12 


< 


1.7 


< 


8.5 X 10^^ 


Q(8) 


8773.22 


< 


0.6 


< 


1.9 


X IO12 


< 


2.9 


< 


8.7 X IO12 










< 


9.9 


X IO12 






< 


3.4 X 10" 



6- value = 1.0 km s ^. 
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TABLE 4 

Column Densities for HD 37903 



Species 


Wavelength 


vlsr 


Nobs 




(A) 


(km s~^) 


(cm~^) 


N{CR) 


4300.313 


10.7 


(9.2 ± 0.1) X 10^^ 




3878.774 




< 4.1 X 10^^ 




3886.409 


11.8 


(2.5 ± 0.7) X 10^^ 




3890.217 




< 9.6 X 10^^ 








(9.2 ± 0.1) X 10^^ 


JV(CH+) 


4232.548 


9.2 


(1.3 ± 0.1) X 1013 




3957.692 


10.4 


(1.2 ± 0.1) X 10" 


Ntot{CE+) 






(1.3 ± 0.1) X 10" 


JV(CN) 


3874.00 


9.3 


(2.3 ± 1.0) X 10" 




3874.61 


10.1 


(5.6 ± 0.9) X 10" 




3875.76 




< 8.8 X 10" 








(7.9 ± 1.5) X 10" 


N{Ca II H) 


3968.468 


-22.9 


(7.4 ± 0.7) X 101° 


7V(Ca II K) 


3933.663 


-25.4 


(1.6 ± 0.1) X 10" 


Af-25(Ca II) 






(1.2 ± 0.1) X 10" 


jV(Ca II H) 


3968.468 


-13.0 


(2.1 ± 0.1) X 10" 


N{Ca II K) 


3933.663 


-14.7 


(2.8 ± 0.1) X 10" 


#-15 (Ca II) 






(2.5 ± 0.1) X 10" 


JV(Ca II H) 


3968.468 


-1.6 


(2.6 ± 0.1) X 10" 


JV(Ca II K) 


3933.663 


-3.3 


(2.6 ± 0.1) X 10" 


iV_2(Ca II) 






(2.6 ± 0.1) X 10" 


JV(Ca II H) 


3968.468 


6.2 


(4.6 ± 0.1) X 10" 


7V(Ca II K) 


3933.663 


5.1 


(4.8 ± 0.1) X 10" 


W+5(Ca II) 






(4.7 ± 0.1) X 10" 


W(Ca II H) 


3968.468 


11.5 


(2.3 ± 0.1) X 10" 


W(Ca II K) 


3933.663 


10.5 


(2.0 ± 0.1) X 10" 


#+10 (Ca II) 






(2.1 ± 0.1) X 10" 


JV(Ca II H) 


3968.468 


16.0 


(8.2 ± 0.5) X 10" 


JV(Ca II K) 


3933.663 


18.1 


(8.9 ± 0.6) X lOi" 


#+18 (Ca II) 






(8.5 ± 0.4) X lOi" 


#to((Ca I) 


4226.728 


8.8 


(4.8 ± 1.5) X 10^ 


#-25 (Na I) 


5895.924 


-20.7 


(2.2 ± 2.7) X 10" 


#-10 (Na I) 




-10.5 


(2.2 ± 0.4) X 10" 


#-5(Na I) 




-4.3 


(2.1 ± 0.4) X 10" 


#+4(Na I) 




3.9 


(2.6 ± 0.5) X 10" 


#+10 (Na I) 




10.0 


(1.1 ± 0.2) X 10" 


#+20 (Na I) 




19.7 


(5.0 ± 0.7) X 10" 


Nt„t{K I) 


7698.974 


10 


1.1 X 10" " 



"Column densities determined from weighted average of detections. 

''6-value = 1.0 km s^^ for CH, CN, Ca I; 2.5 km s^^ for CH+ and Ca II; 2.0 km s^^ for Na I and K I. 
"^Welty; private communication. 
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TABLE 5 
Column Densities for HD 147009 



Species 


Wavelength 


VLSR 






(A) 


(km s~^) 


(cm~^) 


7V(CH) 


4300.313 


4.1 


(1.5 ± 0.1) X 10^^ 




3878.774 




< 3.3 X 10^^ 




3886.409 




< 1.9 X 10^^ 




3890.217 




< 4.4 X 10^^ 


NtotiCB) 






(1.5 ± 0.1) X 1012 


7V(CH+) 


4232.548 


2.8 


(2.4 ± 0.1) X lO^^ 




3957.692 


2.4 


(2.0 ± 0.2) X 10" 


Ntot{C}l+) 






(2.3 ± 0.1) X 10" 


JV(CN) 


3874.61 




< 5.5 X 10" 








< 8.3 X 10" 


JV(Ca II H) 


3968.468 


1.0 


(3.4 ± 0.8) X 10" 


7V(Ca II K) 


3933.663 


1.8 


(3.3 ± 0.2) X 10" 


NtotiCa. II) 






(3.3 ± 0.2) X 10" 


Ntot{Ca I) 


4226.728 




< 3.2 X 10" 


Ntoti^a I) 




1.9 


(1.5 ± 0.3) X 10" 


NtotiK I) 


7698.974 


2.2 


(2.3 ± 0.1) X 10" 



"Column densities determined from weighted average of detections. 

•'fi-value = 1.0 km s^i for CH, CN, Ca I; = 2.5 km s^i for CH+ and Ca II; 1.5 km s^i for Na I and K I. 
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TABLE 6 
lUE ARCHIVE FILES 



Star 


V 


Sp. Type 




E{B - Vj 


Image # 


Notes 


HD 31726 


6.15 


Bl V 


-0.21 


0.05 


SWP 8165; LWR 7098 


1 


HD 37023 


6.70 


BO. 5 Vp 


0.09 


0.37 


SWP 4281; LWR 3782 


2 


HD 37903 


7.82 


B1.5 V 


0.07 


0.35 


SWP 06953, 08058, 8055, 8057, 21293-4; LWR 7025 


2,3,4 


HD 55857 


6.11 


BO. 5 V 


-0.26 


0.02 


SWP 14339; LWR 10953 


5 


HD 64802 


5.49 


B2 V 


-0.19 


0.05 


SWP 14308; LWR 10939 


1 


HD 147009 


8.05 


AO V 


0.26 


0.38 


SWP 22733-4; LWR 17366-7 


2 


C Aql 


2.99 


AO V 


0.01 


0.01 


SWP 8269-70; LWR 7292-4 


3,6 


20 Aquilae 


5.37 


B3 V 


0.09 


0.33 


SWP 08974, 08976, 13865, 36939-40 


4 


HD 200775 


7.20 


B2 Ve 


0.34 


0.45 


SWP 09836-7, 25734, 54089 


4 



^ Unreddened star spectra, average of the two used as comparison for HD 37903 (Fitzpatrick & Massa 1990). 
^Reddened star spectra used for extinction curves (Fitzpatrick & Massa 1990). 
^Witt, Bohlin, & Stetcher 1981. 
"NEWSIPS data used for UV analysis. 

^Unreddened star spectra used as comparison for HD 37023 (Fitzpatrick & Massa 1990). 
^Unreddened star spectra used as comparison for HD 147009. 
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TABLE 7 

UV Results for HD 37903 and HD 200775 











HD 37903 




HD 200775 


Absorption 


Wavelength 


/ 


WA(obs) " 






Wx(ohs) " 


WA(fit) 


Feature 


(A) 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


i^CO (0 - 0) 


1544.448 




215 ± 15 




191 ± 43 






Cl - 0) 


1 50Q 74Q 


0.0291 


43 ± 14 


33 


162 ± 15 


ioO It Z 1 


91 1 


(2-0) 


1477.565 


0.0401 


27 ± 10 


37 




1 no -U '^"7 




(3-0) 


1447.352 


0.0346 


35 ± 10 


34 


26 ± 23 


Q7C -I- 
o / o ± 01 


OQQ Q 


(4-0) 


1419.044 


0.0233 


28 ± 8 


27 


24 ± 10 


111 =t 10 


1 7R 
1 / D 


(5-0) 


1392.525 


0.0145 


24 ± 9 


19 


58 ± 29 






(6-0) 


1367.623 


0.00795 








1 / U Hi oD 


140 


(7-0) 


1344.186 


0.00414 








1 nn -1- 1 n 
luy it ly 


117 
11 / 


(8-0) 


1322.150 


0.00200 








1 1 Q _1_ Qfi 


QO 7 

oy. / 


(9-0) 


1301.403 


0.00095 








QO _|_ 1 c: 


oo 


(10 - 0) 


1281.866 


0.00041 








Q1 _|_ Ofl 
Ol ± ZD 


oZ.O 


(11 - 0) 


1263.433 


0.00018 








fi7 -1- Of\ 
D / ± ZU 


fil A 
01. 4 


i^CO (1 - 0) 


1510.407 


0.0351 








ZDD it oU 


1 
IZu 


(2-0) 


1478.793 


0.0402 








1 QQ -I- Q1 

loo it ol 


1 
IZO 


(3-0) 


1449.145 


0.0347 








1 c: -I- OQ 
IZO It Zo 


1 1 Q 
llo 


(4-0) 


1421.340 


0.0242 








1 1 O _1_ Q7 

IIZ It ol 


1 HQ 
lUo 


(5-0) 


1395.229 


0.0145 








117 -1- 1 R 
11/ =t la 


yo 


C I (J = 0) 


1276.483 


0.00168 




3 


41 ± 13 


1 yi _i_ 1 n 

14 =t lU 


0/1 Q 
Z4.0 




1280.135 


0.0243 


21 ± 6 


20 




90 ± 22 


32.9 




1328.833 


0.0580 


34 ± 7 


29 


51 ± 11 


27 ± 5 


46.2 




1 'ifiO 309 


0.0804 


34 ± 5 


38 


56 ± 14 


D4 ± 14 


oy.o 


(J = 1) 


1276.750 


0.00239 


15 ± 6 


7 


28 ± 9 


no _i_ 1 n 

ZZ =t lU 


00 

zy.z 




1279.056 


0.00202 


11 ± 5 


6 




-1- 1 1 
oZ ±11 


Ofi fi 
Zo.O 




1279.890 


0.0137 


15 ± 4 


18 




CO _|_ 1 Q 

OZ Hz la 


Qf^; Q 

00.0 




1280.597 


0.00684 


14 ± 6 


12 




KQ -i- 01 
Oo It Zl 


QQ 
oo 




1329.085 


0.0193 


45 ± 5 


54 


70 ± 13 


HQ -i- 1 
Do it / 


1 0/1 
1Z4.Z 




1560.682 


0.0603 


72 ± 7 


60 


67 ± 20 


(17 _1_ 1 O 

y / it iz 


lino 

1 iU.Z 


(J = 2) 


1277.723 


0.0145 


9 ± 4 


7 


30 ± 15 




30.7 




1277.954 


0.000967 


5 ± 4 


1 




25 ± 15 


21.2 




1279 229 


00324 








1 Si -1- R 
lo It O 


OQ R 
ZO.O 




1279.498 


0.000903 








QO _L 1 1 
OZ ± 11 


on Q 
zU.o 




1280.847 


0.00505 








Qfl _|_ 1 K 

OD Hi 10 


Of^ A. 
Z0.4 




1329.578 


0.0435 


18 ± 5 


21 


41 ± 12 


(?Q _t_ n 
Do it y 


7t^ Q 
{ O.o 




1561.340 


0.0121 


40 ± 8 


34 




1 ofi -1- 1 n 
izD It ly 


1 1 o c: 

iiy .0 


Blends 


1277.245 


0.0967 


70 ± 7 


62 


78 ± 11 


113 ± 16 


130.9 




1277.513 


0.0242 


41 ± 7 


45 


36 ± 11 


100 ± 19 


96 




1280.333 


0.0147 


18 ± 5 


15 




77 ± 16 


60.2 


S I 


1295.653 


0.087 


10 ± 3 


9 


15 ± 7 


8.8 ± 4.2 


21.9 




1316.543 


0.0381 


12 ± 6 


5 




19.9 ± 6.6 


11.3 




1425.030 


0.141 




12 


18 ± 9 


25.7 ± 5.9 


26.3 




1425.188 


0.0269 








13.8 ± 6.1 


10.8 




1472.971 


0.0218 








57.2 ± 24 


8.0 




1473.994 


0.0803 




7.6 


30 ± 15 


10 ± 5 


20.6 




1474.379 


0.0179 




2.0 


28 ± 14 


11 ± 5 


7.6 




1807.311 


0.096 


17 ± 7 


15 


69 ± 12 


39 ± 5 


33.5 


Ni II 


1317.217 


0.146 


12 ± 3 


14 




28 ± 10 


29.7 




1370.132 


0.144 


13 ± 4 


15 


22 ± 7 


50 ± 9 


31.2 




1454.842 


0.0516 


9 ± 4 


8 




28 ± 9 


27.1 




1709.600 


0.0666 








31 ± 10 


34.9 




1741.549 


0.0776 


15 ± 5 


15 


36 ± 8 


41 ± 7 


36.8 




1751.910 


0.0638 


13 ± 5 


13 


24 ± 9 


31 ± 6 


35.6 



"IVE NEWSIPS data. 
Joseph et al. 1986. 
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TABLE 8 

UV COLUMN DENSITIES FOR HD 37903 AND HD 200775 



Species 


HD 37903 

JV(fit) 
(cm-2) 


HD 200775 
JV(fit) 
(cm^^) 


JV(„((i^CO) 


(0.9 - 1.1) X lO^" 


(1.9 - 2.1) X 10" 


NWCO) 




(1.0 - 1.2) X 10" 


NtotiC I (J = 0)) 


(1.1 ± 0.5) X 10" 


(3.2 ± 1.0) X 10" 


Ntat(C I (J = 1)) 


(1.5 ± 0.5) X lO^" 


(3.6 ± 1.3) X 10" 


Ntot(C I (J = 2)) 


(3.5 ± 1.2) X 10^^ 


(3.2 ± 1.5) X 10" 


Ntot{?> I) 


(6.8 ± 2.4) X 10^2 


(2.6 ± 1.0) X 10^3 


W/„/(Ni II) 


(1.1 ± 0.6) X 1013 


(1.0 ± 0.,3) X lO'" 



TABLE 9 

RATES AND RATE CONSTANTS FOR THE CHEMICAL ANALYSIS 



Reaction Rate/Rate Constant Value Notes 



C+ + CH ^ + H ki 3.0 X 10-1° cm^ g-^ 1 

C2 + CO + C fcj 4.0 X 10-11 cm'' s-i 2,3 

C2 + N -> CN + C h 1.7 X 10-11 cm'' s-i 3,4 

CH + N -> CN + H ki 2.0 x 10-" cm" s-i 3,5 

C+ + NH -> all products h 2.8 x lO-i" cm'' s-i 4 

CN + -j- CO + N ke 1.8 x IQ-" cm" s-i 3,6 

C2 + -s- 2 C GoiC-i) 2.0 x 10-1" g-i 7 

CN + to ^ C + N Go(CN) 1.0 X 10-1° s-i 8 



iPederman & Huntress 1989. 
^Estimate. 

"Rate constant has a temperature dependence of (T/300)°'^. 
"Federman et al. 1994. 
'^Messing et al. 1981. 
"Schmatjko & Wolfrum 1976. 

^Pouilly et al. 1983, corrected for appropriate radiation field; van Dishoeck 1987. 
^Lavendy et al. 1987. 



-38- 



TABLE 10 

RESULTS OF THE CHEMICAL ANALYSIS 



Physical Parameters 



HD 37903 



HD 147009 



HD 200775 



n (cm ^) 
T(K) 

E(B-V) 

Tuv 

7V(CH)„6, (cm-2) 
N{C-2)oi,s (cm-2) 
N{C2)pr,d (cm-2) 
JV(CN)„6, (cm-^) 
iV(CN)^„d (cm-2) 



500 - 1500 
40 

1.0 - 3.0 
0.35" 

2.9'' 

(9.2 ± 0.1) X 10^ 
< 9.9 X 10^^ 

4.0 X 10^^ 
(7.9 ± 1.5) X 10^ 

8.5 X 10" 



400 - 1200 
40 

1.0 - 3.0 
0.38" 

1.7= 

(1.5 ± 0.1) X 10^^ 

< 3.4 X lO^^ 
6.5 X 10^^ 

< 8.3 X 10" 
7.4 X 10" 



300 - 900 
40 

1.0 - 3.0 
0.45" 
4.05'' 
(3.2 ± 0.2) X 10^=''' 
1.2 X lO^^"* 
1.2 X 10^^ 
(1.2 ± 0.3) X lO^^"* 
3.9 X 10^^ 



"For HD 37903 and HD 200775 (Buss et al. 1994) and for HD 147009 (Racine 1968). 
''r„„ = 2A„, but R„ = 4.11 (Buss et al. 1994), somewhat larger than average. 
"^Extinction is similar to the value toward p Oph, r„„ scaled by E{B — V). 
"•Paper I. 
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TABLE 11 

RESULTS OF THE ATOMIC ANALYSIS 



Physical Parameters 



C Per 



HD 37903 



HD 147009 



HD 200775 



n (cm ^) 
G„(Na I) s-i 
G„(K I) 

7V(Na (cm-2) 
N{K (cm-2) 
JV(„j(H) (cm-2) 
(Na I/K 
(Na I/K 

(Na Ijfor/Na Per)pred 
(Na Is4„r/Na Pcr)o6s 
(K Is(ar/K I,; Per)preil 
Per^obs 



1200 ° 
4.5 X 10- 
1.8 X 10- 
(7.3 ± 0.3) X 
(7.3 ± 0.3) X 
(1.6 ± 0.3) X 
67 ' 
100 ± 6 



IQ13 c 

10" = 
10" 



500 » 
1.2 X 10-" 
4.4 X 10-" 
(1.1 ± 0.2) X 10^^ 

1.1 X 10" 
(3.6 ± 0.7) X 10" 
62 

100 ± 18 
0.23 ± 0.06 
0.15 ± 0.03 
0.21 ± 0.06 
0.15 ± 0.01 



400 " 
2.9 X 10-1° 
9.7 X 10-1" 
(1.5 ± 0.3) X 10^^ " 
(2.3 ± 0.1) X 10" " 
(1.9 ± 0.4) X 10" 
56 
65 ± 13 
0.01 ± 0.01 
0.21 ± 0.04 
0.01 ± 0.01 
0.32 ± 0.02 



300 ' 
3.4 X 10 
1.1 X 10 
(5.9 ± 0.8) X 
(1.1 ± 0.1) X 
(2.3 ± 0.3) X 
54 " 
54 ± 9 ' 
0.05 ± 0.01 
0.08 ± 0.01 
0.06 ± 0.01 
0.15 ± 0.02 



11 b 
10 6 
1012 b 
IQlO b 

10" " 



»For 4„ = 1.0. 
'^Paper I. 

"6-value = 1.5 km s'^ for ( Per and HD 147009; 6-value = 2.0 km s'^ for HD 37903. 

Assumes a 20 % error in A'io((H). 
"Includes improvements and corrections to Paper I; see text for details. 
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TABLE 12 
Independent determination of luv/n 



Parameters HD 37903 HD 147009 HD 200775 

(luvin) (cm^) 2.0 x lO"'* " 2.5 x lO"'* " 3.3 x lO"'* ' 

Ntot{^)ab> (cm-2) 3.6 X 10^1 " 1.9 x 10^^ 2.3 x 10^^ ' 

N{R2)„hs (cm-2) 1.0 X 10^1 " 2.5 x 10^° ' 3.2 x 10^° " 
(Iuv/n)H2 (cm^) 0-12 0.09 0.10 



=Prom Tabic 10. 

''Federman ct al. 1997, with corrections; see text for details. 
=Buss ct al. 1994. 
''From Table 11. 
"Vallee 1987. 
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TABLE 13 
UV Results for 20 Aquilae 



Absorption 


Wavelength 


/ 


W\(ohs) " 


Wx(&t) 


Wx ' 


Feature 


(A) 




(mi) 


(mk) 


(mk) 


i^CO (0 - 0) 


1544.448 




150 ± 15 




150 ± 10 


(1 - OJ 


1509.749 


0.0291 


113 i 11 


109 


115 ± 5 


(2-0) 


1477.565 


0.0401 


121 ± 9 


113 


130 ± 5 


(3-0) 


1447.352 


0.0346 


114 ± 9 


107 


120 ± 5 


(4-0) 


1419.044 


0.0233 


106 ± 8 


96 


110 ± 6 


(5-0) 


1392.525 


0.0145 


84 ± 8 


84 


90 ± 8 


(6-0) 


1367.623 


0.00795 


69 ± 8 


70 


75 ± 10 


(7-0) 


1344.186 


0.00414 


59 ± 8 


57 


63 ± 6 


(8-0) 


1322.150 


0.00200 


39 ± 8 


43 


45 ± 10 


(9-0) 


1301.403 


0.00095 


25 ± 6 


30 




(10 - 0) 


1281.866 


0.00041 


16 ± 4 


18 




(11 - 0) 


1263.433 


0.00018 


12 ± 4 


10 




i^CO (1 - 0) 


1510.407 


0.0351 


97 ± 9 




95 ± 6 


(2-0) 


1478.793 


0.0402 


37 ± 4 


38 


31 ± 5 


(3-0) 


1449.145 


0.0347 


25 ± 6 


34 


23 ± 6 


(4-0) 


1421.340 


0.0242 


29 ± 5 


27 


25 ± 8 


(5-0) 


1395.229 


0.0145 


19 ± 8 


20 


12 ± 6 


CI (J = 0) 


1276.483 


0.00168 


19 ± 4 


14 


20 ± 5 




1280.135 


0.0243 


42 ± 6 


43 


45 ± 5 




1328.833 


0.0580 


56 ± 6 


52 


43 ± 6 




1560.309 


0.0804 


63 ± 7 


64 


80 ± 10 


(J = 1) 


1276.750 


0.00239 


7 ± 4 


7 


11 ± 5 




1279.056 


0.00202 


22 ± 5 


6 


23 ± 5 




1279.890 


0.0137 


22 ± 6 


22 


30 ± 5 




1280.597 


0.00684 


15 ± 5 


15 


25 ± 5 




1329.085 


0.0193 


56 ± 6 


61 


70 ± 10 




1560.682 


0.0603 


58 ± 6 


64 


75 ± 5 


(J = 2) 


1277.723 


0.0145 


22 ± 4 


16 


25 ± 5 




1277.954 


0.000967 


7 ± 3 


2 


7 ± 5 




1279.229 


0.00324 


8 ± 3 


5 


12 ± 5 




1279.498 


0.000903 


3 ± 3 


2 


5 ± 4 




1280.847 


0.00505 


16 ± 5 


7 


20 ± 5 




1329.578 


0.0435 


34 ± 5 


41 


40 ± 6 




1561.340 


0.0121 


53 ± 6 


60 


45 ± 5 


Blends 


1277.245 


0.0967 


98 ± 7 


87 


120 ± 10 




1277.513 


0.0242 


68 ± 5 


62 


60 ± 8 




1280.333 


0.0147 


32 ± 6 


27 


30 ± 8 


S I 


1295.653 


0.087 


31 ± 8 


21 


20 ± 10 




1316.543 


0.0381 


31 ± 8 


24 






1425.030 


0.141 


27 ± 6 


24 


20 ± 8 




1425.188 


0.0269 


16 ± 4 


17 






1472.971 


0.0218 


20 ± 7 


16 






1473.994 


0.0803 


23 ± 6 


22 


15 ± 10 




1474.379 


0.0179 


14 ± 6 


15 






1807.311 


0.096 


29 ± 4 


31 


35 ± 10 


Ni II 


1317.217 


0.146 


14 ± 6 


13 


15 ± 10 




1370.132 


0.144 


12 ± 5 


14 






1741.549 


0.0776 


15 ± 6 


14 


20 ± 10 



"lUE NEWSIPS data. 
''Hanson et al. 1992. 
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TABLE 14 

UV COLUMN DENSITIES FOR 20 Aquilae 



Species 




7V(fit) 
(cm-2) 






(4.25 - 4.5) X 10i=> 






(1.8 - 2.0) X IQi" 


NtotiC I (J = 


0)) 


(7.2 ± 2.0) X lO^" 


NtotiC I (J = 


1)) 


(2.4 ± 0.6) X 10^'' 


NtotiC I (J = 


2)) 


(9.5 ± 2.7) X 101== 


NtatiS I) 




(1.1 ± 0.2) X 10^^ 


JV(„,(Ni II) 




(1.0 ± 0.4) X 10^^ 
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TABLE 15 
UV Results for X Per 



Absorption 


Wavelength 


/ 

J 


Wx ( GHRS) 




W\{IUE] 


Feature 


(A) 




(mA) 


(mA) 


(mA) 


C I (J = 0) 


1260.734 


0.0394 


37.1 ± 2.7 


39.1 


44.3 ± 5.5 




1276.483 


0.00168 


28 ± 1.1 


27.4 


29.4 ± 7.2 




1277.245 


0.0967 


42.0 ± 1.6 


40.5 


74.6 ± 6.6" 




1280.135 


0.0243 


38.1 ± 1.5 


39.5 


34.6 ± 6.4 


(J = 1) 


1260.927 


0.0131 


25.8 ± 1.8 


25.3 


46.1 ± 8.1' 




1260.996 


0.00984 


25.4 ± 1.8 


24.3 






1261.123 


0.0164 


27.7 ± 2.0 


26.5 






1276.750 


0.00239 


17.0 ± 0.9 


16.8 


16.8 ± 4.1 




1277.282 


0.0725 


35.7 ± 1.3 


36.6 






1277.513 


0.0242 


26.1 ± 1.0 


27.1 


49.5 ± 6.0^^ 




1279.056 


0.00202 


15.6 ± 0.9 


15.4 


9.6 ± 4.1 




1279.890 


0.0137 


26.0 ± 1.0 


26.7 


20.7 ± 5.1 




1280.404 


0.00438 


21.5 ± 0.9 


21.0 


28.9 ± 5.5 




1280.597 


0.00684 


24.2 ± 1.0 


23.4 


12.9 ± 3.8 


(J = 2) 


1261.426 


0.00984 


11.8 ± 1.2 


14.5 


35.2 ± 10.7'* 




1261.552 


0.0295 


20.1 ± 1.5 


20.0 






1277.5.50 


0.0812 


26.3 ± 1.0 


25.8 






1277.723 


0.0145 


16.0 ± 0.9 


17.0 


16.1 ± 5.5 




1277.954 


0.000967 


2.1 ± 0.6 


2.7 






1279.229 


0.00324 


8.7 ± 0.7 


8.5 






1279.498 


0.000903 


1.9 ± 0.6 


2.1 






1280.333 


0.0147 


17.8 ± 0.9 


17.3 


34.0 ± 4.3" 




1280.847 


0.00505 


10.4 ± 0.8 


11.1 





''Blend with A1277.282. 
''Blend with A1260.996. 
"^Blend with A1277.550. 
■'Blend with A1261.552. 
=Blend with A1280.404. 



TABLE 16 
Ground-based Results for 20 Aquilae 



McD (Moderate Res)" McD (High Res)" McD (High Res)' CTIO (Moderate Res) Other 



Species 


Wavelength 












1>hSK 










(A) 


(mA) 


(km s-i) 


(mA) 


(km s-i) 


(mA) 


(km s-i) 


(mA) 


(km s-i) 


(mA) 


CH 


4300.313 


13.0 ± 0.2 


3.8 


13.3 ± 0.3 


2.8 


14.2 ± 0.7 


2.4 


14.3 ± 0.6 


3.3 


12.6 ± 1.3° 




3878.774 


1.5 ± 0.2 


1.4 


















3886.409 


4.0 ± 0.3 


1.6 










3.1 ± 0.8 


0.9 






3890.217 


3.0 ± 0.3 


2.1 










< 2.3 






CH+ 


4232.548 


2.8 ± 0.2 


4.0 


2.6 ± 0.3 


1.3 


2.8 ± 0.3 


5.0 


2.6 ± 0.7 


0.2 


3.3 ± 0.8'^ 




3957.692 


1.7 ± 0.2 


2.6 










2.9 ± 1.4 


0.2 




CN 


3874.00 


3.2 ± 0.2 


1.3 


3.2 ± 0.3 


2.3 






< 3.7 








3874.61 


11.5 ± 0.2 


0.7 


11.3 ± 0.3 


2.2 






13.4 ± 2.3 


1.1 


10.0 ± 2.0'^ 






1.6 ± 0.2 


0.6 


1.7 ± 0.3 


2.5 






< 3.7 






Ca II 


3968.468 


2.9 ± 0.7 


-15.0 










5.5 ± 3.4 


-20.4 








13.4 ± 0.5 


-8.9 










28.9 ± 4.3 


-8.6 








61.5 ± 0.7 


0.9 










62.3 ± 3.8 


0.7 






3933. dd3 


3.4 ± 0.9 


-18.5 


5.9 ± 1.1 


-14.0 


4.6 ± 0.7 


-15.6 


4.1 ± 0.1 


-18.3 








32.6 ± 0.8 


-7.9 


25.2 ± 0.2 


-7.7 


24.7 ± 0.6 


-8.4 


43.5 ± 2.4 


-9.2 












33.8 ± 0.8 


-0.2 






82.0 ± 2.2 


0.1 








85.2 ± 0.8 


2.1 


37.2 ± 0.8 


2.8 


77.0 ± 1.2 


1.3 
















10.1 ± 1.0 


7.3 












Cal 


4226.728 


0.8 ± 0.1 


-5.4 


0.7 ± 0.1 


-7.3 
















4.2 ± 0.2 


3.5 


3.7 ± 0.2 


1.7 


2.9 ± 0.3 


2.9 


2.1 ± 0.5 


3.5 


4.0 ± 0.2'^ 


Na I 


5895.924 


89.0 ± 35.8 


-8.0 






78.0 ± 1.0 


-7.7 












210.6 ± 43.6 


-1-2.4 






186.0 ± 2.0 


-1-2.3 










5889.951 


94.8 ± 28.5 


-8.1 




















244.4 ± 41.0 


-1-2.4 
















K I 


7698.974 


6.2 ± 1.3 


-8.9 






4.5 ± 0.6 


-7.8 












111.3 ± 1.4 


-1-2.6 






110.0 ± 1.0 











''2dcoude. 

''6 foot camera. 

'^Hanson et al. 1992. 
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TABLE 17 

C2 DATA TOWARD 20 AqUILAE 



Line 


X 


Wi 


» LSR 


N 


Notes 




(k) 


I'm A1 


(km s^^) 


(Wl j 




R(0) 


8757.69 


2.0 ± 1.0 


7.4 


(3.0 ± 1.5) X 10^^ 


1 






2.0 ± 0.2 




(3.0 ± 0.3) X 10^^ 


2 


N(0) 








(3.0 ± 0.3) X 10^^ 


3 


P(2) 


8766.03 


< 1.1 




< 1.6 X lO^'' 


2 


Q(2) 


8761.19 


3.9 ± 0.5 


7.2 


(1.2 ± 0.2) X 10^^ 


1 






7.3 ± 0.5 




(2.3 ± 0.2) X 10" 


2 


R(2) 


8753.95 


3.4 ± 0.7 


6.8 


(1.3 ± 0.4) X 10 


1 






5.0 ± 0.4 




(1.9 ± 0.2) X 10" 


2 


N{2) 








(1.8 ± 0.1) X 10" 


3 


P(4) 


8773.43 


1.8 ± 0.3 




(1.6 ± 0.3) X 10" 


2 


Q(4) 


8763.75 


3.3 ± 0.4 


7.8 


(9.9 ± 1.3) X IO12 


1 






5.7 ± 0.5 




(1.8 ± 0.2) X 10" 


2 


R(4) 


8751.68 


3.1 ± 1.4 


5.2 


(1.3 ± 0.3) X 10" 


1 






3.2 ± 0.2 




(1.5 ± 0.1) X 10" 


2 


iV(4) 








(1.4 ± 0.1) X 10" 


3 


P(6) 


8782.31 


3.1 ± 0.6 




(2.4 ± 0.5) X 10" 


2 


Q(6) 


8767.76 


4.6 ± 1.6 


8.5 


(1.4 ± 0.5) X 10" 


1 






4.2 ± 0.4 




(1.3 ± 0.1) X 10" 


2 


R(6) 


8750.85 


< 4.95 




< 2.5 X 10" 


1 






3.0 ± 0.3 




(1.5 ± 0.2) X 10" 


2 


N(6) 








(1.4 ± 0.1) X 10" 


3 


Q(8) 


8773.22 


< 1.8 




< 5.4 X IO12 


1 






1.1 ± 0.3 




(3.3 ± 0.9) X IO12 


2 


N(8) 








(3.3 ± 0.9) X 10" 


3 


Ntot{C2) 








(5.2 ± 1.6) X 10" 


4 



^6-value = 1.0 km s^^. 

^Federman et al. 1994. 

^Weighted average of individual detections. 

''Sum of the weighted averages of the 0,2,4,6,& 8 rotational Hnes. 
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TABLE 18 
Column Densities for 20 Aquilae 



Species 


Wavelength 


VLSR 


Nobs 




(4) 


(km s^i) 




JV(CH) 


4300.313 


3.0 


(2.0 ± 0.1) X 10^^ 




3878.774 


1.4 


(1.1 ± 0.1) X 10^^ 




3886.409 


1.6 


(9.4 ± 0.7) X 10^^ 




3890.217 


2.1 


(1.1 ± 0.1) X 10" 


JV(„,(CH) 






(2.0 ± 0.1) X 10" 


JV(CH+) 


4232.548 


2.3 


(3.2 ± 0.2) X 10^^ 




3957.692 


1.7 


(3.8 ± 0.5) X 10^^ 








(3.3 ± 0.2) X 10^^ 




0074.UU 


1.8 


/I 1 _L n 1 ^ 1 nl2 
(1.1 ± U.IJ X iU 




3874.61 


1.3 


(3.1 ± 0.1) X 10^2 




3875.76 


2.1 


(1.1 ± 0.1) X 10^2 








(4.2 ± 0.1) X 10 


JV(Ca II H) 


3968.468 


-17.7 


(2.3 ± 0.7) X 10" 


JV(Ca II K) 


3933.663 


-16.6 


(4.8 ± 0.1) X 10" 


JV_i7(Ca II) 






(4.8 ± 0.1) X 10" 


JV(Ca II H) 


3968.468 


-8.8 


(3.4 ± 0.1) X 10" 


JV(Ca II K) 


3933.663 


-8.3 


(3.6 ± 0.1) X 10" 


JV_8(Ca ii) 






(3.5 ± 0.1) X 10 


JV(Ca II H) 


3968.468 


0.8 


(2.4 ± 0.1) X 10" 


JV(Ca II K) 


3933.663 


-0.2 


(5.0 ± 0.2) X 10" 


7Vo(Ca II) 






(2.9 ± 0.1) X 10" 


7Vi(Ca II K) 


3933.663 


1.2 


(2.2 ± 0.1) X IO12 


JV3(Ca II K) 




2.8 


(5.6 ± 0.2) X 10" 


JV7(Ca II K) 




7.3 


(1.2 ± 0.1) X 10" 


JV7(Ca I) 


4226.728 


-7.3 


(2.8 ± 0.3) X 10" 


iV3(Ca I) 




2.9 


(1.4 ± 0.1) X 101° 


7V(Na I Di) 


5895.924 


-8.0 


(1.0 ± 0.1) X 10" 


JV(Na I D2) 


5889.951 


-8.1 


(1.8 ± 4.0) X 10" 


iV_8(Na I) 






(1.0 ± 0.1) X 10" 


JV(Na I Di) 


5895.924 


2.4 


(1.4 ± 0.1) X 10" 


JV(Na I D2) 


5889.951 


2.4 


(7.3 ± 3.0) X 10" 


iV3(Na I) 






(1.3 ± 0.1) X 10^^ 


Af-9(K I) 


7698.974 


-8.9 


(2.8 ± 0.4) X 10" 


#3(K I) 




2.6 


(1.4 ± 0.1) X 10" 



^Column densities determined from weighted average of detections. 

''6-value = 1.0 km s^i for CH, CN, Ca I, Na I and K I; 2.5 km s^i for CH+ and Ca II. 
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Fig. 1. — Representative interstellar atomic spectra toward HD 37903 are shown here. Top 
panel — Ca II K, solid line is high resolution 2dcoude data; dashed line is moderate resolution 
2dcoude data. In all figures the fiux is normalized to unity. There is a nice correspondence 
between the two difTerent resolutions in velocity and line shape. Center panel — moderate 
resolution 2dcoude Ca II H data. Bottom panel — moderate resolution 2dcoude data; dashed 
line is Na Di and sohd line is Na D2. The slight velocity differences are attributed to an 
error in wavelength solution. 
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Fig. 2. — Representative interstellar molecular spectra toward HD 37903 are shown. Top 
panel — CH, solid line is high resolution and dashed line is moderate resolution 2dcoude 
data. Center panel — CH+, same as top panel. Bottom panel — CN, same as top panel. 
There is a general agreement in velocity and hne shape; the slight velocity differences are 
attributed to an error in wavelength solution. 
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Fig. 3. — Moderate resolution 2dcoude molecular spectra toward HD 37903 arc shown. Top 
panel — CH, solid line is A — X transition and dashed line is B — X transition. Center panel 
— CH"*", solid line is ^4 — X (0—0) transition and dashed line is A — X (1—0) transition. 
Bottom panel — CN R(0) transition is shown. 
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Fig. 4. — Representative interstellar atomic spectra toward HD 147009 are shown. Top panel 
— Ca II K, dark solid line is high resolution 2dcoude data; dark dashed line is moderate 
resolution 2dcoude data; and blue dashed line is moderate resolution CTIO data. Center 
panel — Ca II H, dark dashed line is moderate resolution 2dcoude data; and blue dashed 
line is moderate resolution CTIO data. Apparent minimum in the CTIO data is noise, not 
an additional component. Bottom panel — shows spectra of several atomic species taken 
with different instruments at different resolutions. There is excellent agreement among all 
observations. Solid line is K I and dark long dashed line is Na Di high resolution 6 foot 
camera data; dark dashed line is moderate resolution 2dcoude Na Di data; and pink dashed 
line is moderate resolution 2dcoude Na D2 data. Use the left scale for K I and the right scale 
for Na I D lines. 



-SI- 



LOS 
1.00 
0.95 - 
0.90 - 
0.85 - 



1.05 

>r 1.00 
w 

g 0.95 

I 0.90 h 
> 

'a 0.85 
0.80 
0.75 
1.1 

1.0 
0.9 
0.8 
0.7 - 
0.6 - 



I I I I I I I I 

HD 147009 




CH ?i4300 

I I I I I I I I I I 



McDonald 
Cerro Tololo 



I I I I I I I I I I I I I I I 




CH* ?i4232 

I I I I I I I I I I 



I I I I I I I I I I I I I I I 



V 



Kl ?i7699 



1 



1.25 
1.00 

- 0.75 

- 0.50 
0.25 
0.00 



-30 



-20 



-10 



10 



20 



30 



Vlsr (km s" ) 



Fig. 5. — Representative interstellar spectra toward HD 147009 are shown. Top panel — 
CH, dark solid line is high resolution 2dcoude data; dark dashed line is moderate resolution 
2dcoude data; and blue dashed line is moderate resolution CTIO data. Center panel — CH+, 
solid line is high resolution 2dcoude data; dark dashed line is moderate resolution 2dcoude; 
and blue dashed line is moderate resolution CTIO data. The slight velocity differences are 
attributed to an error in wavelength solution. Bottom panel — shows spectra of K I, solid 
line (left scale), and Na Di, dashed line (right scale), taken with the high resolution 6 foot 
camera. 
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Fig. 6. — Interstellar molecular spectra toward HD 147009 are shown; dark lines are 2dcoude 
data and blue lines are CTIO data. Top panel — CH, the solid lines are the A — X transition 
and the dashed line is the B—X transition. Center panel — CH"*", the solid lines are the A—X 
(0—0) transition and the dashed lines are the A — X (1—0) transition. The slight velocity 
differences are attributed to an error in wavelength solution in these moderate resolution 
data. Bottom panel — shows the high resolution 2dcoude 3874 A region of the spectrum 
where CN resides; no detection is evident. 
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Fig. 7.— a.) Representative NEWSIPS lUE ^^CO spectra toward HD 37903 are shown. 
The data are represented by the filled circles. Our best fit to the data (solid line) and the 
data— fit (dashed line, offset to 1.22) are also shown. 
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was not fit. Data— fit offset to 1.30. 
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Fig. 8. — b.) Same as figure 8a, but spectra of ^^CO are shown. The e — X (5—0) band of 
^^CO is clearly seen in the spectrum as well. 
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Fig. 9. — Shown are the derived extinction curves for the reddened hnes of sights. The 
dashed curve is the average interstellar extinction curve represented by ( Per scaled to the 
E{B — V) of the reddened star. Extinction curves for ( Per, 20 Aquilae, and HD 200775 
were obtained from Savage & Mathis 1979, Savage et al. 1985, and Walker et al. 1980, 
respectively. 
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Fig. 10. — Results of C I excitation for HD 37903 and HD 200775 are shown. Density is for 
collisional partners, see text. The dashed hnes — A^(J=1)/A^(J=0) ratio and the dashed- 
dotted hnes — A^(J=2)/A^(J=0) ratio based on ± 1 a uncertainty in column densities. For 
HD 37903 (top panel), the arrow indicates upper range for A^(J=1)/A^(J=0) ratio is a lower 
limit. Results for = 3 are indistinguishable for luv = 1- For HD 200775 (lower panel), 
upper limit for iV(J=2)/A^(J=0) ratio yields much larger density. Blue scale shows results 
for 3-fold enhancement in radiation field where distinguishable. Result from C2 excitation 
is also shown. 
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Fig. 11. — b.) Same as Figure 11a, but spectra of ^^CO are shown. 
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Fig. 12. — Representative atomic spectra toward 20 Aquilae are shown. Top panel — 
Ca II K, dark solid line is high resolution 6 foot camera data; dark long-dashed line is 
high resolution 2dcoude data; dark dashed line is moderate resolution 2dcoude data; and the 
blue dashed line is moderate resolution CTIO data. Center panel — Ca II H, dark dashed 
line is moderate resolution 2dcoude; blue dashed line is moderate resolution CTIO data. 
Bottom panel — shows spectra of several atomic species taken with different instruments 
at different resolutions. There is excellent agreement among all observations. Solid line is 
K I and dark long-dashed line is Na Di high resolution 6 foot camera data; saturated dark 
dashed line is moderate resolution 2dcoude Na Di data while weaker dark dashed line is K I; 
and pink dashed line is moderate resolution 2dcoude Na D2 data. 



- 64 - 



1.05 

1.00 

0.95 

0.90 

0.85 

0.80 

0.75 
1.02 

>. 

4-« 

■« 1.00 

0) 

4-« 
C 

§! 0.98 

'.5 
« 

0.96 



1.00 
0.98 
0.96 
0.94 
0.92 

-30 -20 -10 10 20 30 

Vlsr (km s"") 

Fig. 13. — Representative spectra toward 20 Aquilae are shown. Top panel — CH, dark solid 
line is high resolution 2dcoude data; dark long-dashed line shows the high resolution 6 foot 
camera data; dark dashed line is moderate resolution 2dcoude data; and blue dashed line is 
moderate resolution CTIO data. Center panel — CH"'", solid line is high resolution 2dcoude 
data; dark long-dashed line shows the high resolution 6 foot camera data; dark dashed line 
is moderate resolution 2dcoude data; and blue dashed line is moderate resolution CTIO 
data. Velocity offset is due to error in wavelength solution. Bottom panel — shows spectra 
of several atomic species taken with different instruments at different resolutions. There is 
excellent agreement among all observations. Solid pink line is Ca I high resolution 2dcoude 
data; pink and dark long dashed line shows Ca I and Na Di high resolution 6 foot camera 
data; and dark solid line is K I high resolution 6 foot camera data. Use the right scale for 
Na I and K I and the left scale for Ca I. 
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Fig. 14. — Moderate resolution molecular spectra toward 20 Aquilae are shown; dark lines 
are 2dcoue data and blue lines are CTIO data. Top panel — CH, solid lines are A — X 
transition and dashed lines are B — X transition. Center panel — CH+, solid lines are A — X 
(0—0) transition and dashed lines are A — X (1—0) transition. Bottom panel — CN R(0) 
high resolution 2dcoude data is shown. 



